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CE) 株式 会 社 TREND-PRO / 漫画 绘制 
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一 提 到 微分 方程 ， 大 家 总 感觉 挺 难 。 事 实 上 ， 这 门 课 的 确 挺 难 懂 。 老 实说 ,我 当 
时 听课 的 时 候 ， 也 总 是 搞 不 明白 。 正 如 本 书 中 的 野山 大 地 所 说 ， 虽 然 会 计算 ， 但 是 完 
全 不 知道 为 什么 要 那么 算 。 尽 管 记忆 了 各 式 各 样 的 解 题 模式 和 公式 ， 练 习题 也 能 解 得 
出 来 ， 但 很 多 人 对 微分 方程 总 有 一 种 云 山 雾 单 的 感觉 。 

一 般 大 家 都 认为 求解 微分 方程 是 一 件 不 太 容易 的 事 ， 但 本 书 中 不 会 出 现 解 不 出 
的 微分 方程 。 知 道 解法 的 话 ， 谁 都 应 该 能 解 得 出 来 。 数 学 家 们 已 经 用 各 种 解法 和 公式 
为 我 们 开辟 了 许多 条 道路 ， 沿 着 这 些 道路 前 进 ， 谁 都 可 以 得 到 解答 。 只 是 ， 艰 难 跋涉 
于 不 熟悉 的 数学 世界 时 ， 我 们 往往 会 被 眼前 的 数学 运算 难 住 。 如 果 仅仅 圈 于 自己 的 小 
圈子 ， 就 无 法 通 览 全 局 ， 最 后 只 会 迷失 在 数学 的 世界 里 。 抬 头 环顾 一 下 四 周 ， 也 许 你 
会 发 现 另 一 番 开阔 美丽 的 景色 。 

这 本 书 是 一 本 循循善诱 的 微分 方程 方面 的 指导 书 ， 它 会 带 你 畅游 在 微分 方程 的 世 
界 中 。 与 一 般 的 教科 书 有 所 不 同 ， 本 书 没有 网 罗 微分 方程 这 门 学 问 的 所 有 分 支 ， 也 没 
有 刻意 追求 数学 的 严密 性 和 普遍 性 。 总 之 ， 在 沿 着 既定 路 线 艰难 摸索 的 同时 ， 也 请 你 
尽情 地 享受 沿途 的 美景 吧 。 就 像 我 们 头顶 上 的 天 空 一 样 ， 独 自 飞翔 于 微分 方程 的 世界 
里 也 是 有 很 多 乐趣 的 。 人 类 虽然 没有 翅膀 ， 但 是 已 经 能 够 制作 翅膀 在 空中 生 翔 了 。 有 
了 微分 方程 这 个 翅膀 ， 我 们 也 一 定 能 够 在 数学 世界 里 自由 地 畅 翔 。 如 果 大 家 将 本 书 作 
为 跳板 ， 也 能 在 微分 方程 的 天 空中 展翅 高 飞 的 话 ， 我 会 荣幸 之 至 。 

最 后 ， 本 书 能 够 面世 多 亏 了 欧姆 社 的 全 体 员工 、 用 “数学 之 神 ” 的 有 趣 创意 为 本 
书 创作 脚本 的 SWP 先 生 、 成 功 地 将 抽象 数学 世界 表现 为 生动 形象 图 画 的 漫画 家 徐 子 ， 
在 此 向 大 家 表示 衷心 的 感谢 。 本 书 的 面世 是 我 们 团队 合作 的 结果 。 
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一 一 不 只 是 摇摆 运动 





网 啡 的 冷却 
2. 火箭 的 飞行 . 






























































2 \ 
$ | 好 无 聊 啊 ! 
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数 宫 神社 的 守护 神 
数学 女神 ( 年 龄 不 详 ) 




















































让 我 借用 一 下 你 的 
有 ! 身体 吧 ， 


es oo 



























































































































































































































































































































































































































































































































































erates 
神社 里 没有 神仙 ae 
恐怕 不 行 吧 ? zy 

但 是 ， 又 没有 来 

SHE, AME 

WR! 



















































































































































































































































































































































































































































































































































































































































































ATRAI | 。 最 近 一 点 利润 都 没有 ， 
MARETE? 



















































我 有 义务 守护 这 个 神社 。 


身 为 学 问 之 神 ， 还 请 您 
好 好 工作 才 是 ! 



























































































































































































































































































































































































































































































































































































































好 ! 
实现 ! 
绝对 让 它 实现 ! 















































































































































那么 ， 请 您 满足 下 一 个 
SHAH BBC! 





























Mp FAI LMA, 
NY 但 也 不 错 。 | | 
| 


0 
o 





l [f | 大 学 = 年 级 学 生 
野山 大 地 





| FF, RAR! 


困惑 的 年 
BA! 










































































































































































那 也 太 便 宜 他 了 吧 












































































































































您 为 什么 要 诗 这样 pE > 
A BBR? i 


AS ae 
x. 即便 会 计算 ， 也 不 明 
自 自己 在 做 什么 。 


所 以 ， 你 想 更 好 更 深 
入 地 去 理解 它 …… 
































干什么 去 呀 ? 




























































































































































































BH, BK, GRRL 
AURA, HUT 


还 是 先 学 习 吧 


你 …… 不 太 了 解 微分 方 
程 ， 倒 是 对 甜点 比较 熟 












































































































































ot | 



















用 它 可 以 模拟 地 操纵 
飞机 。 


住 的 现实 世界 的 现象 。 


那么 ， 飞 行 模拟 装置 是 如何 
感知 飞机 的 飞行 的 呢 ? 


Sh 
虽然 还 不 是 
PPTA WwW 
有 和 和 


大 概 














































































































微分 方程 是 一 种 工具 ， 能 够 
在 数学 世界 里 表现 你 们 所 居 























































































这 可 是 微分 方程 的 
功劳 啊 ! 





首先 用 数学 的 抽象 概 
念 表达 现实 世界 的 飞 


通过 求解 微分 方程 ， 
能 够 得 到 函数 。 

















顺便 说 一 下 ， 将 现象 
模型 化 之 后 得 到 的 数 
学 模型 是 包含 微分 的 












































能 够 完成 模拟 实验 ， 
当然 也 能 用 于 预测 。 





em y 
EANA.) 














现实 世界 REER 























BIRARE 


s 


















































将 大 气 的 运动 模型 化 ， 得 到 机 
分 方程。 通过 求解 这 个 殷 分 方 bo 
可， 就 能 名 所 测 出 大 气 的 活动 





























除 此 之 外 ， 微 分 方程 还 被 广泛 应 微分 方程 也 可 用 于 预测 今后 全 
用 于 灾害 规模 、 人 口 变化 、 广 告 4 < 球 变星 的 情况 以 及 全 球 变 暖 可 
效果 和 商品 销售 额 等 各 种 各 样 的 能 引起 的 结果 。 

领域 。 



































































































































































































































































































































= - ie 原来 如 此 啊 ! 























































































































举 个 实际 的 例子 ， 
怎么 样 ? 
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说 什么 呢 ! 



























































































































































如 何 得 知 那 架 滑翔 机 
























































































































































































































































































































































@ 将 现实 中 滑翔 机 的 
运动 模型 化 

















加 利用 这 个 函数 ， 我 们 可 以 得 
知 盖 翔 机 从 过 去 到 未 来 的 任意 
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我 们 来 了 解 一 下 关于 £ 
物体 运动 的 一 些 定律 TREAT! 
和 法 则 。 Ses 

a 


m (mass) = 物体 的 质量 
a (acceleration) = 物体 的 加 妾 度 
F (force) = 作用 力 


位 置 对 时 间 的 微分 形式 表示 
加 速度 ， 


那么 ， 也 就 是 说 ， 能 够 用 
微分 方程 表达 物体 的 运动 
方程 。 





顺便 说 一 下 ， 含 徽 分 的 A 
方程 就 是 微分 方程 。 这 i 这 样 我 们 或 完成 了 滑翔 


机 运动 的 简单 模型 化 。 


通过 微分 方程 ， 我 们 就 能 得 知 
滑翔 机 从 过 去 到 未 来 任意 时 刻 




































































































































































ELBA OM RE 


提炼 出 这 个 本 质 的 模型 就 什 
得 感动 …… 








ak 
那样 啊 …… 
















































































































































































































































































对 于 数学 女神 来 说 
微分 方程 一 写 出 来 


就 明白 所 有 的 事 
情 了 。 


麻烦 您 ， 能 帮 我 求解 
这 个 微分 方程 吗 ? 


实 世界 的 壁 全 ， 大 地 
也 跨 不 过 数学 世界 的 
EAN. 


{ 数学 大神 难以 越过 现 





那么 ， 
dr ` 
a ea ee 








这 次 ， 我 们 求 一 下 满足 
这 个 运动 方程 



























































求解 微分 方程 ， 要 得 


到 的 不 是 数值 ， 而 是 


不 能 像 求 方程 的 解 那样 ， 
代入 数值 进行 求解 ! 





除法 一 乘法 


为 求解 微分 方程 ， 
必须 进行 积分 。 
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滑翔 机 以 恒定 的 过 度 在 空中 
飞行 的 情况 下 ， 机 身受 到 的 
作用 力 是 保持 平衡 的 。 








ENBER GLAD ASR 
A, Eh BERTARA GLY 
分 力 与 浮力 分 别 保 持平 衡 。 


s : 为 求解 方程 ， 设 0 与 
也 就 是 说 ， 所 受 合力 : TARK, HER 
为 替 ， 因 此 运动 方程 
PH AFAR, 





在 表示 时 间 的 变量 t 中 ， 


代入 想 得 知 滑翔 机 的 过 去 位 置 
或 者 未 来 的 某 个 时 刻 。 


就 可 以 得 到 与 这 个 时 刻 相 
对 应 的 滑翔 机 的 位 置 。 
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a(t)=ut+2, 





























MER) 


























































































































现 家 的 说 明 





























如 果 得 到 的 结果 发 生 严 
重 偏离 ， 就 再 一 次 返回 
到 模型 化 这 个 步骤， 考 
处 更 加 精密 的 模型 。 





微分 方程 模型 














徽 分 方程 模型 





























采用 何等 精度 的 模型 取决 
于 起 从 多 远 的 距离 晓 望 肖 
Faw, 


若 想 知道 滑翔 机 周围 空气 的 
流动 状态 ， 应 该 考虑 更 加 复 
杂 和 精密 的 模型 。 








BAP, REHM 
望 滑翔 机 整体 的 运动 时 ， 可 
以 忽略 其 形状 ， 当 做 一 个 点 
来 处 理 。 


















































因为 在 这 种 情况 下 ， 不 ss 
能 忽视 机 身 的 形状 
模型 化 是 经 验 ! 为 充分 观察 现实 世界 中 发 生 
的 现象， 
oe 
pees eT Ti 


ABH RH FOR, 













































































































































































































































































对 于 典型 的 现象 ， 可 以 考虑 
各 种 各 样 的 模型 ， 因 此 ， 和 要 
学 会 使 用 这 些 模型 。 


类 似 的 问题 都 需 
要 会 解 才 行 。 





从 今 以 后 ， 你 要 好 好 学 习 
贯通 现实 世界 与 数学 世界 
的 本 事 。 


\ ee 知 …… 知 道 了 ! 








好 好 努力 呀 ! 
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€ 
el / 
E> 
fs 
保佑 我 在 这 次 模 NG 
拟 考 试 中 考 出 好 eye y Sits Bde + 


idi j ERE, RRA 
keee 














= E 学习 去 ! M Z ŽIA! Si Rodis 






































我 不 是 学 这 个 参拜 者 又 不 多 ， 怎 么 
这 样 对 他 们 呢 ， 


虽说 他 们 看 不 到 ， 但 是 也 
HELM ER! 




















人 们 把 烦恼 什么 的 
部 依赖 在 我 身上 。 


但 是 听取 这 些 不 就 
是 数学 女神 您 的 主 
要 工作 吗 ? 





我 早 就 到 了 。 
War 
ce 类 
BD 
在 里 面 帮忙 打扫 人 
来 着 。 m 


















































我 想 知道 ， 为 什么 能 
微分 方程 进行 模型 


用 
化 


i 
分 了 吗 ? | 


Z 


看 来 ， 学 习 典 型 微分 方程 
的 模型 之 前 ， 要 先 钻研 重 
要 基础 知识 才 行 啊 













这 些 概念 在 前 几 天 的 讲解 
中 都 已 经 用 到 了 ， 在 模型 
HH EAH HR, 





设 速度 为 "， 温 度 为 T， 


但 是 ， 将 现象 模型 
化 时 ， 经 常用 文字 
RAR, AAS. 
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OEE TEET ET 
面孔 了 ! 


研究 一 下 与 现实 世界 相对 应 


的 数学 世界 的 语言 。 









语 
MNI : 
YNE) P thafoy ap i, 
iy RHEE. 这 样 的 ay 
py $i A. A 
变量 是 将 任意 的 No ° 


a 





数 用 文字 表达 的 

















































































































WR HbA S 
po TI 比如 说 ， 写 成 Jz) 的 这 


个 式 子 表示 变量 V 是 变量 








Z 的 值 随意 变化 时 ， 与 
此 相对 应 ，y 的 值 也 发 
生变 化 。 



































即使 是 相同 的 变量 ， 
8 和 z 的 情况 也 是 不 对 
等 的 。 


























































































































| 真 像 现在 的 我 和 Cray) 
因此 ，z 岂 做 自 变量 ， 你 的 关系 呀 ! 
4 叫做 因 变 量 。 PE 
yuo f 
» op ¢ x 











| TN 
自 变量 与 因 变量 是 
主 从 关系 。 






























































































































































































































































ee 像 不 像 水 森 RL, ARF, HA MHS 
和 数学 女神 的 关 动 模型 化 时 必要 的 变量 ! 坐标 用 来 表示 位 置 关系 ， 
AF? 经 常用 >、%、2z 等 符号 来 
表示 。 
假设 滑翔 机 沿 着 工 轴 前 
it, 
对 不 起 


— 
j; o FANAKA G 
op iz 




























































































































































































假设 现在 滑翔 机 以 25m/s 的 迷 


度 飞行 ， 与 时 钟 上 显示 的 数值 
相对 应 的 i 时 间 ) 轴 上 
那么 可 以 用 滑翔 机 的 位 置 的 数 表示 时 间 。 
随时 间 的 变化 表示 滑翔 机 
的 运动 。 








位 置 随时 间 变 化 还 是 时 间 随 
位 置 变化 呢 ? 






那么 ……, 考 虑 滑翔 机 运动 模型 的 | 


函数 时 ， 必 须 确定 Z 和 t 当 中 ， 哪 
个 是 自 变 量 ， 哪 个 是 因 变 量 。 







哪个 随 着 哪个 变 


了 能 够 清楚 地 表现 主 从 
时 间 t 是 自 变量 ， a M aT Rta 
位 移 x 是 因 变 量 .。 达 习惯 来 表现 。 


还 是 选择 位 置 对 于 
时 间 的 关系 吧 。 
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这 样 一 未 …'… 利 用 画 数 
符号 ， 可 以 用 z=f(t) 





这 是 位 置 对 于 时 间 的 函 
数 ， 因 此 用 7 替换 /就 
能 用 z=z(t) 


AARE- TAMA | kant 
运动 。 25m/s 的 速度 在 
飞行 ， 





Z 的 单位 为 来 ， 
t 的 单位 是 秒 ， 
因此 可 以 表示 
0 成 z(t)=25t。 



































































































































那么 用 曲线 表示 滑翔 
机 的 位 置 随时 间 的 变 








就 变 成 这 样 一 条 直线 。 


该 曲线 表达 的 是 滑翔 机 
的 位 置 z 对 于 时 间 t 的 函 
数 z(t)。 


作曲 线 的 话 ， 你 们 
EERE BS 
数 的 法 则 。 


从 过 去 到 未 来 的 时 间 
里 ， 每 过 去 1 秒 钟 ， 滑 
期 机 的 位 置 就 会 远离 
25 来 。 
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另外 ， 设 直线 的 斜率 a 是 Az 与 At 的 比值 ， 即 Az/At。 其 
Ra, 中 ，At 为 直线 在 横 坐 标 上 的 变化 
量 ，Az 为 直线 与 At 相对 应 的 纵 坐 


标 上 的 变化 重 ， 


是 希腊 文字 中 的 德 


了 愤 重 起 见 在 这 里 
Re! 为 了 慎重 起 见 在 这 里 声明 


一 下 ， 由 于 A 不 是 变量 ， 
所 以 不 能 约 分 。 






























































应 的 是 这 里 的 速度 ， Å 当 直 线 的 针 率 更 陕 的 情况 下 ， 
pn RARE RES ABR RAT ZSMS. 



























































































































































是 25m/s 对 不 对 ? 
xim] 
Ke X(t)=25t/m 
125 
100 
aS 
50 
25- 
ol t $ t + + t CS) 








那么 ， 如何 从 曲线 上 得 到 迷 
度 的 具体 数值 呢 ? 


是 这 样 的 ， 保 留 单位 的 


情况 下 ， 计 算 位 移 对 于 
时 间 的 变化 率 Az/At。 


若 时 间 t 的 单位 为 s( 秒 )， 位 务 z 的 
单位 为 m( 米 )， 那 么 速度 的 单位 就 
是 m/s( 米 每 秒 )。 


未 每 小 时 
nee 


顺便 说 一 下 ， 曲 线 的 截 距 j 是 数学 世 
界 的 语言 ， 表 示 曲 线 与 纵 坐 标 相 交 
的 数值 。 


在 这 个 例子 中 ， 与 b 对 应 的 现实 
世界 中 的 量 是 初始 位 移 。 


ROG, ARABIA 
期 机 所 处 的 位 置 。 





若 时 间 的 单位 为 h( 小 时 )， 位 移 的 
单位 为 km( 千 来 )， 那 么 速度 的 单 
位 就 是 km/h( 千 米 每 小 时 )。 





X(t)=25t+b 
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著 不 用 具体 的 数值 ， 为 了 路 入 数学 世 
界 ， 要 合 弄 一 些 
个 性 才 行 啊 ! 


而 是 用 文字 像 Z=at+b 


LAAT PR RS, 
则 更 具 普 遍 性 。 





不 只 是 这 些 ! 滑翔 机 的 
tR, +è, afti 






































REAM LETH 
比如 说 ， 如 果 用 25m/s 表 +, 
FARMER, RAR 
达 一 个 飞行 的 状态 。 2 
数学 世界 中 表现 相同 
的 事物 者 能 够 用 同样 
的 方法 进行 处 理 








但 用 wv 表示 过 度 的 话 ， 而 
能 够 表示 任意 的 迷 度 。 



















































































这 也 是 数学 世界 的 优点 
eI 





























进入 数学 世界 ， 





就 没有 必要 介意 
这 些 变 量 原来 表 
示 的 是 什么 。 


在 现实 世界 中 解释 数学 世界 得 
出 的 结果 时 ， 必 须要 把 个 性 重 
新 找 回来 才 行 啊 ! 





alae Ti 起 起 这 里 的 变量 表示 


的 是 什么 。 








数学 世界 中 有 一 些 是 
BEA HH HR, È 
们 的 名 字 分 别 被 附 在 
书 的 后 面 。 








感觉 有 点 明白 模型 化 








要 能 够 闪 认 其 中 一 
些 常用 的 函数 及 其 
曲线 。 
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图 指数 函数 





| 时 而 连续 递增 、 时 而 连续 递减 





设 o 为 某 一 常数 ， 则 表示 成 
y(z) =a" 


的 函数 称 为 指数 函数 。 y=0.5" y=0.257 so 


指数 函数 的 曲线 如 右 图 所 示 。 


指数 函数 中 ，a 为 任意 常数 ， 在 处 理 微分 “ 
方程 时 非常 重要 的 指数 函数 是 以 


lim ( + 1J- e= 2.718281828459045235--- 


定义 的 常数 e 的 指数 函数 。 

y(a)=e 
这 里 出 现 的 符号 e 表 示 的 数 叫 纳 皮尔 数 。 指 数 函 数 的 基 
本 性 质 有 


s, z+y Q z- z) zI a 本 
daaa”, = (EPan, (§) = 


等 等 。 


Ws OF 05108 


© (x) Be 





BAZ% 
EETL 








指数 函数 的 反 函 数 称 为 对 数 函数 。 我 们 把 Y y=log, 
T= 
改写 成 y 等 于 的 式 子 
y(x) = log x 
这 就 是 对 数 函数 。 
y= log, 952 


对 数 函 数 的 曲线 如 右 图 所 示 。 


常数 a 称 为 底数 ， 底 数 为 10 的 对 数 log1oz 是 常用 对 
数 ， 底 数 为 e 的 对 数 log.z 称 作 自然 对 数 。 自 然 对 数 在 处 理 
微分 方程 方面 尤其 重要 。 为 了 表示 自然 对 数 ， 可 以 简化 
为 
y(a) = Ine Vz)-Inz 的 曲线 
对 数 函 数 的 基本 性 质 有 
log,cy = log x + logy, log, 3 =log,2 -log,h, log,z"= ylog,z, log, = re 由 


等 等 。 
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EZAN 
| 时 而 递增 时 而 递减 摇摆 变化 的 函数 








如 下 图 所 示 ， 以 直角 三 角形 的 各 个 边 长 与 角 的 关系 
i g b 
sing =—, COST = — 
a a 


定义 的 函数 统称 为 三 角 函 数 。 由 于 这 些 函 数 是 以 2 为 周期 的 周期 函数 ， 所 以 常常 用 
它们 来 描述 周期 性 的 现象 。 根 据 勾 股 定理 ， 有 


coss + sin’a = 1 


成 立 !。 设 p=cosz，q=sinz， 则 P+ 2=1， 表 示 与 pg 平面 圆周 上 的 点 相对 应 。 





© 三 角 函数 的 曲线 一 一 左边 的 是 sinz、 右 边 的 是 cosz 





1 利用 接 下 来 要 出 现 的 欧 拉 公式 ， 也 可 以 得 到 ， 
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ENHAN 
| S=mamiems memes 








乍 一 看 没有 任何 关系 的 指数 函数 与 三 角 函 数 ， 通 过 欧 拉 公 式 ? 
Ct"=cosz ising 


引入 虚数 形成 难以 分 离 的 关系 。 用 欧 拉 公式 


ny ei _ eet 


cose, sing = i 
可 以 得 到 上 述 关系 式 。 双 曲线 函数 的 定义 与 此 非常 相似 。 双 曲线 函数 是 定义 为 
Cret ipn Or 
7 , sinhr= 7 
的 函数 ( 右 图 ) 。 根 据 定义 ， 
eret Y (e-e Ý 
content (SE (e) 


„etei _ ete? 
4 4 


coshz = 





=1 


> 双 曲线 函数 的 曲线 一 一 左边 为 sinhz， 


成 立 。 设 p=coshz，g=sinhz， 则 p?-q?=1 ， 与 pg 平 pe 


面 双 曲 线 上 的 点 相对 应 。 这 就 是 双 曲 线 函 数 这 个 

名 字 的 缘由 。 由 于 双 曲线 函数 被 称 作 hyperbola， 因 此 cosh 被 称 为 余弦 双 曲 线 函数 ， 
sinh 被 称 为 正弦 双 曲线 函数 。 虽 然 ， 名 字 与 三 角 函 数 颇 为 相似 ， 但 是 函数 的 性 质 方面 
基本 没有 相似 的 地 方 ， 所 以 还 是 把 它们 当做 完全 不 同 的 东西 看 待 吧 。 





2 将 指数 函数 与 三 角 汕 数 分 别 按 级 数 展开 ， 可 以 得 出 如 下 推论 。 当 z=x 的 时 候 ， 变 成 e*+1=0，e、x、1、i、0 结 合 在 一 起 , 得 到 了 一 
个 神秘 的 式 子 。 无 理 数 e 的 ix 次 方 加 1 得 9， 实 在 是 令 人 惊讶 。 
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但 是 ， 实 际 上 也 能 
减速 和 加 速 的 吧 ? 








刚才 考虑 了 滑翔 机 以 恒 
定 的 速度 飞行 的 情况 。 





一 直 以 恒定 的 速度 持 
续 飞行 的 滑翔 机 之 类 
的 情况 当然 不 存在 ! 













TE: 


R, EREET 
也 不 能 着 陆 ! 










所 以 ， 为 了 使 模型 更 加 接近 
现实 世界 ， 我 们 来 考虑 一 下 
运动 状态 变化 的 情况 吧 ! 






假设 飞行 路 线 为 直线 ， 
过度 是 变化 的 。 









UR MARE, RAH EBL 
为 纵 坐 标 作 图 可 以 得 到 曲线 。 














这 种 速度 变化 的 情况 RRA TABS Bee, 
TF, Rae A RA 就 是 说 可 以 用 变化 率 AT/At 表 
FARR AHF? 示 














速度 恒定 的 情况 下 ， 图 表 为 
直线 因此 不 用 考虑 变化 率 ， 


SRRELHHR, BARE 
成 曲线 ， 因 此 当时 间 的 变化 
量 较 大 时 ， 与 此 相对 应 的 变 
化 率 也 会 发 生变 化 。 
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这 时 ， 我 们 试 着 
ils as ee AN 
} | 


变化 率 则 表示 曲线 上 那 
个 点 的 切线 的 斜率 。 


时 间 的 变化 量 At 变 成 0 的 话 ， 
Az/At 中 的 分 母 就 变 成 0， 所 
以 就 无 法 计算 了 。 










不 能 把 At 变 
RO, MAR 
Ath oreo 


像 这 种 ， 使 变量 无 限 
趋 近 于 某 个 数 的 数学 











Al 一 0 


(无限 趋 拓 于 o 的 情况 ) 












不 管 怎 么 说 ， 这 里 能 
够 求 出 的 变化 率 只 有 
没有 理由 就 是 aaa 
最 好 的 理由 ! 





这 样 就 能 够 连接 现实 世 oo 
界 和 数学 世界 了 。 


`g 




















现在 的 情况 是 ， 求 时 间 


















































































































































MAERA TOW 
情况 下 ， 时 间 的 变化 量 
与 其 对 应 的 位 移 变化 量 
的 比值 Az/At 的 极限 。 
Ame —T Rt 可 以 这 样 来 表现 。 
处 理 ! 
除法 计算 前 面 的 
lime 那么 ， 位 置 的 变化 量 Az 
是 表示 从 t 时 刻 的 位 置物 位 置 随时 间 的 变化 
动 到 t+At 时 刻 的 位 置 的 fyr = 
、 li BBE. 公式 (2) a j ( 
-am Az = x(t + At) - x(t) # 
At 一 0 
AAT “AR” TN 
(limit) 的 符号 。 可 以 这 样 表示 。 


用 上 述 减 式 赫 换 位 置 的 变 ARRE JOR: 


化 量 Az 的 话 ， 则 时 间 间 a(t + At) - x(t) 
也 At 无 限 趋 近 于 0 hn = 


这 个 式 子 叫 做 与 时 刻 相 
对 应 的 位 置 古 数 Z(t) 的 -lim Zt AD = 2) 
导数 。 At 
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© 


我 们 习惯 性 地 称 这 种 
数学 运算 为 …… 


Z(t 对 t 求 微分 。 




























































































dt doh? 









































































































































































































































如 果 表 示 为 时 间 的 函数 ， 则 较 
为 麻烦 ， 


所 以 有 时 只 用 ZK 表 示 ， 
d Ab 如 dz/dt。 


dt 


A 























































































































































































































狂想 岗位 音 点 区 示 ， 2 
也 可 以 写成 …… 
四 四 名 
RAGE r 
读 成 “dxdt”。 
表示 函数 zt) 
的 导数 。 
有 
不 是 分 数 呀 ? 由 于 d 不 是 变量 ， 
所 以 不 能 进行 约 
ah 
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接着 ， 我 们 来 看 看 
位 移 对 时 间 的 变化 
PAz/Atho 





从 图 上 也 可 以 看 出 ，dz/dt 表 
wt A HRR, 


而 AZ/At 表 示 的 是 在 At 的 时 间 
间 硼 里 的 平均 速度 。 








与 导数 dz/dt 相 对 应 
的 量 表 示 的 是 位 移 除 
以 时 间 ， 





ML, dradtR 5th aia 
对 应 的 位 移 7 的 导 孝 也 就 
AR, 


因此 速度 就 变 成 了 时 间 的 
函数 VCt)。 











函数 (tb) 叫 做 函数 
a(t) HH Ke, 


















RRS, AA 
与 横 坐 标 平行 。 





RRS WMA Eta 











BARHK, HUBS 
做 曲线 。 TV a 





横 坐 标 为 时 间 t， 纵 
AAR RR hA 
Wife 















t 

















达 度 与 时 间 戌 正 
HH, BRAS 
FRRHER, 

































HRA} F HR 也 可 以 将 上 式 再 次 
对 时 间 进行 微分 。 














iy 
是 将 酒 数 Z(t) 对 时 间 g i 
rene t aml 还 能 再 次 微分 吗 ? 
AAA 
一 g a ZF 
at 3 

















要 注意 ， 这 里 的 2 虽然 
看 起 来 像 二 次 方 ， 但 
并 不 表示 或 者 1 的 二 
KF 





的 次 数 叫 做 导 函 数 
的 阶 数 。 
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© 稍微 偏离 了 滑翔 机 的 运动 ， 开 始 研究 起 微分 计算 的 相关 内 容 了 。 


Q xm 
(人 水木 真是 有 干劲 呵 
QO) wm see ERM. DOL RINE, AUR OR Fe 











数 ft) 的 微分 。 
首先 是 常数 的 情况 
ft) 
根据 微分 的 定义 得 出 
t+At)- ft) _, 1-1 
t f= in £ n £0 im 1 Ao, 
结果 为 0。 
其 次 ， 与 成 正比 的 情况 下 
fit=t 
同样 得 到 
了/-lim KEHAD- Ab) lim CHAP -t lim ttAt-t 1 
dt amo At ao At ao At 
等 于 t 的 平方 的 情况 下 
有 -让 
得 到 
Z io- -lim flt+At)- fit) _ lim +- im (+ 2tAt + At) -È 
At te at ato At 
dim AAEN = lim(2t + At) = 2t 
amo At a 


那么 ， 对 t 求 一 次 导数 得 到 2t。 
不管 怎么 说 能 够 看 到 法 则 之 类 的 东西 了 吧 ? 
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© tb 许 一 列 出来 的 话 ， 会 更 容易 看 伸 。 


d 
J0-1 一 下 和 -0 
d 
人 各 -一 二 人 -1 
d 
sf — = 
At) >  f)=2t 


啊 ! 总 算是 看 出 来 了 …… 


这 里 ， 我 们 使 用 了 数学 里 非常 有 用 的 一 般 化 技巧 。 设 "为 变量 ， 而 不 是 具体 的 数 ， 
考虑 一 下 f 与 t 的 n 次 方 成 正比 的 情况 

fee 
这 时 ， 情 况 会 变 得 稍微 复杂 一 点 ， 可 化 为 


SD 


n+ 的 -他 
dt Tosin At 


pag 
ES (t+At) -t 

At At 

Tim CANETA + + ntAt™ t+ At") —t" 
are) At 


lim nt™ At +--+ ntAt™l+ At" 
ro At 


ant 


最 后 变 成 的 n-1 次 方 的 函数 nt"!， 不 是 吗 3? 
fit)=t" > 2 fit)=nt™ 


©) 这 就 是 所 谓 的 微分 公式 


A ponm , 
dt 


3 式 子 中 出 现 的 省 略 号 表示 的 各 个 项 都 是 存在 的 。 由 于 n 是 变量 ， 无 法 确定 项 数 ， 虽 然 不 能 全 部 写 下 来 ， 但 这 样 写 还 是 可 以 
让 读者 看 慌 ， 所 以 就 这 样 写 了 。 


公式 是 数学 世界 的 便利 工具 。 由 于 有 了 公式 ， 就 没有 必要 和 逐一 从 头 开始 计算 
了 。 话 虽 这 么 说 ， 但 还 是 建议 大 家 试 着 自己 导 一 遍 公式 。 即 使 从 未 独立 导出 过 
一 个 公式 ， 沿 着 前 人 成 功 的 路 线 ， 谁 都 可 以 导出 这 些 公式 。 

谁 都 可 以 ? 

是 的 。 这 个 很 重要 。 计 算 公式 的 步骤 看 起 来 虽然 很 复杂 ， 但 是 所 做 的 事情 既 出 


平 意料 又 简单 。 教 科 书 中 会 告诉 你 几乎 所 有 函数 的 微分 公式 ， 因 此 务必 自己 亲 
手 导出 这 些 公式 为 好 ! 


举 几 个 微分 公式 的 例子 


© @@ © 


A sint = cost 

dt 

d 

— cost = -sint 

dt (2.1) 


fee et 
dt 


而且， 微分 的 基本 性 质 如 下 所 示 。 
d d d 
(oD + Balt) = a5 AO + Br A 


d d d 
yl) = wt) FO) + fa) 


d d d (22) 
qO- ay MO GO 


dz 1 
‘aps sae (y= fiz), z= f(z)) 


dy 
® e gw 
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利用 微分 ， 从 表示 位 移 的 函数 D 
TENA RRR BAU), ‘ 


反 过 来 起 想 从 表示 速度 的 
函数 U(t) 得 到 表示 位 移 的 
六 数 Z(t) 的 方法 吧 ! 


在 加 度 恒定 的 情况 下 ， 从 时 间 
ESNE Ft i F A Atty t; 
里 移动 的 距离 | 


首先 考虑 一 下 以 恒定 
RAAT MTA 
期 机 的 运动 吧 ! 


ix 
=< 














移动 距离 | 正好 是 与 这 个 
长 方形 的 面积 相等 。 


著 能 将 这 个 关系 一 般 化 ， 


什么 样 束 度 的 曲线 都 可 以 
通过 求 面积 得 到 它 的 移动 
距离 ! 
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可 以 推测 是 底 边 长 度 夹 
以 高 度 除 以 2。 


滑翔 机 的 速度 与 时 间 成 
正比 时 ， 


V(t)=at 会 是 什么 
HR? 


这 样 也 不 能 变 成 一 般 化 
的 情况 。 


为 了 求 出 导数 ， 把 时 间 
间隔 At 当做 线索 考虑 。 
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v 
做 法 是 用 过度 恒定 
时 的 方法 ， 求 出 近 
似 值 。 
Bo < 
1 EVALUA FUAT FNAL FUSA 
就 这 样 ， 假 设 固 定 的 
宽度 变 得 足够 细 ， 就 
UW 会 逐渐 接近 v(t)=at。 

















t 
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T =at + WAt + what 
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这 样 做 ， 





把 萝卜 切 成 n 个 长 条 





大 地 真是 能 
vam, ATAR 
这 样 … hal a 
vt 3 
` 
E Di 
g Lav At+v,At+ +v,At= Dv At 
É 


在 字母 表 中 表示 成 5， 


据说 是 求 和 的 意思 。 


RAH, 就 是 对 这 个 
WAt 的 j 按 顺序 从 1 到 n 
得 到 的 所 有 项 的 求 和 。 





为 了 减少 误差 ， 邻 长 
条 形 的 根 数 接近 于 无 
HK! 





暂时 先 把 萝卜 放 
一 边 吧 ! 


积分 符号 了 是 将 字母 














| 胡 中 的 S 上 下 拉 伸 得 


通过 积分 表达 式 ， 就 能 够 
从 表示 过度 的 函数 v(t) 得 
到 表示 位 移 的 函数 z(t) 的 





意思 是 将 v(t) 对 t 从 到 


tj 进行 积分 。 < 


顺便 说 一 下 ，t 称 为 积分 变量 ， 
十 叫 积分 下 限 ，tj 叫 积分 上 限 。 





© 肖 获 机 在 时 刻 的 位 移 坐 标 z(tj) 是 滑 关机 在 时 刻 的 位 移 坐 标 z(b 和 从 4 时 刻 到 时 
刻 的 时 间 段 移动 的 距离 /之 和 ， 则 
a(t,)=a(t,)+1 
若 已 知 速度 的 函数 vfb)， 移 动 的 距离 [等 于 v(b) 对 时 间 ;进行 积分 


x(t =z(6)+ ou 


t 


可 以 表示 成 上 式 。 这 里 ， 为 了 使 静止 ( 表示 特定 的 时 刻 ) 的 变量 动 起 来 (不 只 是 
在 特定 的 时 刻 ) ， 可 以 表示 到 任何 一 个 时 刻 # 的 运动 ， 这 样 一 来 ， 移 动 距离 是 时 
间 的 函数 ， 变 成 


w= ‘wnat (23) 
t 
Be, tnh AA ERER E R a(t) 
a(t)=a(t) + | “uae (24) 
设 位 移 z(t) 满 足 z(t)=0， 以 位 移 为 纵 坐标 、 时 间 为 横 坐 标 ， 作 图 则 如 下 所 示 。 


z 


令 时 间 随 位 置 的 变化 


© 那么 ， 当 得 知 速度 随时 间 变 化 的 函数 w(t) 时， 能 够 求 得 从 t 时 刻 到 任意 时 刻 t 的 时 间 
段 的 移动 距离 ((t) 的 积分 (2.3 ) 也 是 时 间 的 函数 。 


4 或 者 说 任意 的 时 刻 。 
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on 加 


当 积分 上 限 移动 的 时 候 所 求 的 面积 也 会 发 生变 化 
争 求 得 从 t 时 刻 到 任意 时 刻 ! 的 时 间 段 内 移动 距离 !(t) 的 积分 


®© 这 里 ， 将 函数 !(b) 对 t 进 行 微分 。 根 据 微分 定义 可 以 写成 
I(t + At) - I(t) 

At 
FASS RF RUE), W 


d ; 
qO 


Joa -fioa [oa 
à ERT rater 


At 0 At 
它 的 极限 趋 近 于 v(t)。 也 就 是 说 ， 变 成 





4 1) =lim (2.5) 
LO- 
可 是 由 于 函数 tt) 是 式 (23) 
1) = [oat 
积分 得 到 的 函数 再 进行 微分 就 回 到 原 函数 
谷 。 导电 记 相 对 应 的 位 置 z(t) 可 以 表示 为 


a(t)=a(t) + Jroa 
再 一 次 将 表示 任意 时 刻 的 变量 t 设 为 t ， 两 边 都 减 去 z(t)， 左 右边 进行 替换 ， 则 


[par=zt -zt 
ty 


SEARS at 有 这 种 想法 的 你 或 许 对 数学 精通 到 已 经 没有 读 这 本 书 的 必要 了 ， 或 许 就 是 赁 死记 硬 背 来 
对 待 数学 的 。 


这 个 式 子 表示 的 是 ， 滑 翔 机 在 t 时刻 的 位 置 z(t/) 减 去 t 时 刻 的 位 置 z(t) 就 可 以 
得 到 t, 时 刻 到 t/ 时 刻 滑翔 机 移动 的 距离 。 这 是 当然 的 。 但 是 ， 移 动 的 距离 是 v-t 图 
的 面积 。 一 方面 ， 对 位 移 进行 微分 即 可 得 到 速度 这 个 量 。 图 的 面积 和 微分 即 可 得 
到 速度 的 位 移 这 个 量 之 间 理 应 是 没有 任何 关系 的 。 但 是 可 以 用 等 号 将 他 们 连接 起 
来 ,好像 这 里 有 很 重要 的 关系 呀 。 
{ro = fit) (26) 
这 里 的 F(t) 称 为 At) 的 原 函数 。 函 数 Flt) 和 f(t) 的 关系 是 ， 对 函数 Flt) 进 行 微分 
就 可 得 到 导 函 数 ft)， 对 函数 f(t) 进 行 积分 即 可 得 到 原 函 数 Ft)。 就 像 乘 法 和 除法 的 
关系 一 样 ， 微 分 和 积分 是 彼此 的 逆 运 算 。 
‘eh san 
Co ks oo | 
maw 
MERPRIS ERAO, Name FT BL BNA a RRR) 
ORO SRAE AEN 








刚才 的 速度 与 位 置 的 关系 式 
ou-ae ) -z(t) 


用 ft) 和 ft) 的 原 函 数 F(t) 分 别 蔡 换 v(t) 和 z(t)， 用 a 和 6b 分别 替 换 t 和 t,， 则 


ou =F()-F(a) . 


这 就 是 微 积分 的 基本 定理 。 


这 个 式 子 要 表示 的 是 ， 可 以 根据 微分 的 逆 运 算 ( 右边 ) 求 得 函数 的 图 形 所 包围 的 
面积 (左边 ) 。 
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到 此 为 止 说 明 的 是 关于 
徽 分 和 积分 的 微 积分 的 


基本 定理 。 q 
| y 


| ou-ro-zaol 
| 


这 就 是 为 什么 我 希望 你 们 
掌握 机 积 分 的 基本 定理 的 
RA. 









$ 





微分 方程 能 够 预知 未 来 的 
关键 就 在 这 里 。 






不 要 再 摆 出 一 副 完全 
不 知 所 云 的 表情 。 









这 个 定理 可 以 说 
是 人 类 的 财富 ! 



























































左边 求 面积 的 积分 是 着 眼 于 全 
局 的 观点 上 成 立 的 ， 与 此 相对 
应 的 是 六 


右边 的 原 函 数 是 着 限于 图 
片上 某 一 点 的 局 部 观点 成 
的 。 
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大 局 的 观点 ， 就 比 
如 说 为 了 看 某 条 河 
坐 着 滑翔 机 从 上 面 


这 样 也 可 以 吗 ? 大 
局 观点 和 局 部 观点 
在 现实 世界 中 是 完 
全 不 同 的 观点 呀 。 








与 此 相对 应 的 是 ， 局 
部 的 观点 是 从 河 边 观 
察 眼前 的 河流 。 就 是 
这 样 的 意象 。 








也 就 是 说 ， 现 实 世 
界 中 从 河 边 如 何 细 
致 地 观察 河流 ， 


也 看 不 到 河流 的 
ah. 


数学 世界 里 积分 的 大 局 观点 和 向 Q 
分 的 局 部 观点 是 联系 在 一 起 的 。 A 








he c 这 就 是 微 积 分 的 基本 定理 告诉 
` 我 们 的 ! 
















































































sari, an \_P 因此 ， 在 数学 世界 里 ， 从 
学 世界 里 只 知道 S| HAHEAM MERE 
Tht, N TALS RoR AA AL 


就 算 不 能 从 河 边 观察 
到 什么 ， 也 能 保证 看 
ATR SH, 




































































































































































那么 ,从 下 次 开始 学 习 典 
型 的 微分 方程 模型 。 


四 i Åz x 
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微 积分 的 基本 定理 的 意思 是 ， 由 函数 的 曲线 形成 的 图 形 面 积 等 于 原 函 数 中 带 入 积分 
上 限 的 值 减 去 带 入 积分 下 限 的 值 。 用 


LFO- FQ) F(a) 

符号 表示 Fb)Fla)， 则 化 为 
f 

| Kodit =PO 

也 就 是 说 ， 只 要 知道 了 原 函数 ， 求 函数 图 形 的 面积 时 就 不 用 一 一 求 极限 ， 而 是 通过 
原 函数 的 减法 就 可 以 求 得 5。 这 真是 方便 啊 ! 

将 函数 jg) 的 原 函 数 用 如 下 符号 表示 ， 则 

roa 
这 叫做 不 定 积分 !。 与 此 相对 应 ， 求 得 由 函数 曲线 形成 的 图 形 面积 的 积分 为 


b 
Í ftdt 


这 叫做 定 积分 。 就 是 说 ， 微 积分 的 基本 定理 将 看 似 没有 关系 的 不 定 积分 和 定 积分 联 
系 在 了 一 起 。 
那么 ,假设 Ft) 为 f() 的 原 函 数 ，C 为 任意 常数 ， 根 据 微 分 的 性 质 可 以 得 出 


d d 
P+ 0) 5 FO+0= A) 

F(t)+C 也 变 成 了 f(s) 的 原 函数 。 也 就 是 说 ， 不 定 积分 变 成 
|roar= re +o 


常数 C 无 法 确定 。 这 里 的 任意 常数 C 叫 做 积分 常数 。 由 于 积分 常数 无 法 确定 ， 所 以 叫 
做 不 定 积分 。 用 曲线 来 表现 的 话 ， 由 不 定 积分 得 到 的 原 函数 的 曲线 平移 之 后 肯定 能 够 全 
部 重合 。 


6 通过 求 极限 计算 积分 在 现实 中 是 很 难 的 。 
7 注意 这 个 名 称 哦 。 
8 这 里 不 是 指 式 子 的 形式 ， 而 是 指 其 意义 。 
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由 于 积分 是 微分 的 逆 运 算 ， 因 此 将 导 函 数 的 公式 从 相反 方向 读 出 即 可 得 到 不 定 积分 
的 公式 。 


Ding [ae tc (nz0) 
dt n 

 sint =cost 一 J costat=sint+ 0 
L cost=-sint 一 |intar=-eost+ 0 
dy 

—e'=e' = Jeldt=e'+C 

dt 


如 何 具体 地 解释 微分 方程 是 下 一 章 要 讲 的 内 容 ， 下 面 大 致 让 大 家 看 一 下 微分 方程 的 
分 类 ， 然 后 结束 本 章 的 内 容 。 就 像 以 前 看 到 的 那样 ， 微 分 方程 中 只 要 引入 导数 ， 就 会 变 
成 微分 方程 ， 如 果 不 考虑 分 类 的 方法 ， 无 疑 会 杂乱 无 章 、 毫 无 秩序 。 对 微分 方程 进行 分 
类 时 ， 需 要 掌握 根据 自 变量 的 个 数 、 导 函数 的 阶 数 和 次 数 、 线 性 和 非 线性 、 常 系数 和 变 
系数 、 奇 次 和 非 齐 次 等 微分 方程 的 分 类 方法 。 





多 微分 方程 的 分 类 


首先 ， 我 们 从 自 变量 的 个 数 开始 吧 。 一 个 自 变 量 对 应 一 个 因 变 量 的 微分 方程 叫做 常 
微分 方程 。 有 几 个 自 变量 的 情况 下 ， 就 要 引入 偏 微 分 进行 求 导 。 含 偏 微分 的 方程 叫做 偏 
微分 方程 。 本 书 只 讨论 常 微分 方程 。 到 此 为 止 ， 我们 所 见 到 的 微分 方程 也 全 部 都 是 常 微 
分 方程 。 由 于 不 会 出 现 偏 微分 ， 请 把 本 书 中 提 到 的 微分 方程 均 视 为 常 微分 方程 。 

其 次 ,方程 的 阶 数 和 次 数 。 微 分 方程 中 所 含 导 函 数 的 最 大 阶 数 称 为 微分 方程 的 阶 。 


把 这 个 导 函 数 的 次 数 称 为 微分 方程 的 次 数 。 导 函数 的 次 数 表示 的 是 导 函 数 的 寡 数 
( 乘 几 次 ) 。 也 就 是 说 ， 


至 + hz=kha ”+ 一 阶 一 次 微分 方程 


由 于 含 一 阶 导 函数 ， 它 的 次 数 为 一 次 ， 因 此 被 分 类 为 一 阶 一 次 微分 方程 。 另 外 ,在 
moe +y dz + kr-0 o 二 阶 -次 微分 方程 


中 ， 虽 然 导 函数 既 有 一 阶 导 函 数 又 有 二 阶 导 函数 ， 但 是 根据 最 大 阶 数 才 是 微分 方程 阶 数 

的 规则 ， 它 就 是 二 阶 ， 这 个 导 函 数 的 次 数 为 一 次 ， 因 此 就 被 分 类 为 二 阶 一 次 微分 方程 。 
再 次 ， 就 是 线性 与 非 线性 的 分 类 。 线 性 微分 方程 就 是 因 变 量 和 它 的 导 函 数 的 次 数 为 

一 次 的 微分 方程 。 若 含有 不 是 一 次 的 项 ， 就 是 非 线性 微分 方程 。 以 刚才 的 微分 方程 为 例 


hem ka 一 线性 微分 方程 


mt +udz +kz=-0 — 线性 微分 方程 


由 于 因 变量 和 它 的 导 函 数 的 次 数 是 一 次 ， 因 此 都 属于 线性 微分 方程 。 
随后 ， 是 常 系数 和 变 系数 的 分 类 。 线 性 微分 方程 的 系数 全 部 为 常数 时 称 为 常 系数 微 
分 方程 ， 若 系数 是 变量 ， 就 称 为 变 系数 微分 方程 。 比 如 说 


虽然 是 变 系数 微分 方程 ,但 
KERAN N 
ae t at 
却 是 常 系数 微分 方程 。 
最 后 ， 是 齐 次 和 非 齐 次 的 分 类 。 线 性 微分 方程 的 情况 下 ， 与 因 变量 没有 关系 的 常数 


为 0 时 ， 称 为 齐 次 微分 方程 ， 不 为 0 时 叫做 非 齐 次 微分 方程 。 比 如 说 


就 是 非 齐 次 微分 方程 。 
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dz 


mv 


dé? 
为 齐 次 微分 方程 。 


对 微分 方程 进行 分 类 时 ， 要 综合 使 用 上 述 这 些 分 类 方法 。 比 如 说 


dr 


moe gy pend 


dt? 


REM -UERSATCREE ADE. KAKA FATA OR, KERE, AL 
根据 具体 情况 进行 省 略 。 这 一 点 在 处 理 具体 微分 方程 时 要 注意 。 

本 书 分 别 在 第 三 章 讨论 可 分 离 变量 一 阶 齐 次 微分 方程 ， 在 第 四 章 讨论 一 阶 非 齐 次 线 
性 微分 方程 ， 在 第 五 章 讨论 二 阶 线性 微分 方程 。 

















本 书 要 处 理 的 微分 方程 实例 如 下 。 
阶 数 种 类 计算 公式 说 明 章节 
aP 描述 北海 道 放生 存 数字 = 
SERE ae 变化 的 微分 方程 第 三 
一 阶 = 
| me 
| 
dr de 考虑 弹力 和 阻力 的 振动 
齐 次 线性 a a E PE a 体系 的 运动 方程 第 五 章 
二 阶 
非 齐 次 线 dr ,dz 包含 外 力 的 振动 体系 的 
性 me A | 第 五 章 
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这 家 伙 又 来 数 宫 
HRR 





保佑 我 的 店铺 
ll WN sexe: 











a 1 ZA wits 
Se LH J| 得 到 惩罚 的 ! 

















部 说 了 不 是 我 的 
Lh! . 





如 果 能 把 这 个 完成 
的 话 ， 我 也 想 实现 za SAREE DF 
一 下 大 家 的 愿望 ， 


但 是 我 也 不 是 万 能 
的 呀 ! 





S 9 
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神仙 能 够 更 加 具 休 地 
帮助 人 们 实现 他 们 的 
BB. 


我 也 只 是 在 做 自己 
份 内 的 工作 而 已 


偷懒 是 因为 那儿 
的 女 亚 实在 是 太 
烦人 了 




































































仙 的 辛苦 啊 ! 


， ao 
SA = 
p FA 
+ oft j3 \ L serai We Y 
rt 走 AN L4 
7 O 
























































oT 
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今天 ， 我 起 试 着 解释 一 下 
递增 和 递减 的 现象 。 
cd 





Í ama csuroee B 
K fj 


\ | \\\ HY / 





我 们 用 可 分 离 变量 微分 方程 Wit MM 


a ies S 





st ax 真 想 去 啊 ! 








还 是 先 学 习 吧 。 


这 些 东 西 吗 ? 
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车 站 前 面 开 了 一 
家 北海 道 的 零售 
EA, METH 
虚心 。 


那里 卖 的 软 冰 洪 
淋 真 是 绝 了 ! 
































NG 





其 他 的 呢 ? 



































北海 道 斋 剧 增 部 
好 像 要 成 为 一 个 


据说 ， 北 海道 庆 增 加 是 
因为 捕食 北海 道 斋 的 北 A 





BERK TORE 
ÀT. 


f 
l 
Val 
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若 真是 那样 的 话 




















© 
g 
G \ 
ARAMA IE— A a 
| / RAE IE 
Q 型 吧 ! 
\ \\ | I I] 77 
好 N 现实 世界 | 数学 世界 7 Z 
N a Zo 
WA A 


























































































































































































































wR 首先 还 是 要 模型 化 。 


ES HH 
量 的 增加 你 认为 应 该 建立 怎样 
的 模型 呢 ? 





还 不 知道 到 底 怎么 
Woo 





BAF ae 
FARA Re 


新 生 许 的 数量 碱 去 死去 
的 庭 的 数量 就 能 够 得 到 
增加 量 。 


PTEE ET Si 
BREAST, CARTE 
MREP— RR, 活 下 来 的 小 许 的 数量 


并 且 死 亡 的 北海 道 庞 
的 数量 也 与 北海 道 旗 
的 总 数 成 正比 。 


首先 从 简单 的 开始 
考虑 ! 
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我 们 考虑 北海 道 鹿 的 增 
加 量 与 北海 道 订 的 数量 


10 头 北海 道 庞 的 群 里 ， 死 亡 
1 头 ， 新 生 5 头 的 话 。 


100 头 北海 道 鹿 的 群 的 话 
应 该 是 每 年 死亡 10 头 ， 新 生 
50 头 。 


ARTAR, eH 
BR Ayo 


mw 


比例 常数 


u 











顺便 说 一 下 ， 确 定 这 个 函数 


关系 的 变量 叫做 参数 ， 也 叫 
这 里 ， 设 时 间 为 Et， 北 海道 4 Rhee. 


的 生存 数量 为 Ft)， 











E2 
Be se 



































HLT AR A RE AE 明白 了 。 
增加 率 dP(t)/dt 与 北海 道 麻生 存 数 量 k o 
Rb) 成 正比 的 模型 。 wee 9 














那么 ， 式 子 就 变 
成 这 样 。 


























描述 北海 道 放生 序数 顺便 问 一 句 ， 比 例 常 数 /该 
和 量变 化 的 微分 方程 ! ELAR? 
$ A 
~ 
nn 4 
t 1 





调查 一 下 现实 世界 的 情 
况 才 能 知道 啊 ! et = 
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如 果 已 知 局 部 的 状况 ， 那 么 能 
够 从 局 部 的 情况 得 知 全 局 的 情 
况 吧 ? 





~ M 


已 知 某 个 时 刻 的 北海 道 庞 的 生存 数量 | 
量 P 和 决定 增加 方法 的 参数 //， 可 
以 了 解 到 什么 呢 ? 

















这 就 是 为 什么 能 够 预测 到 
一 年 甚至 十 年 后 北海 道 许 
生存 数量 的 原因 。 





用 积分 求 函数 ! 
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现在 ， 我 们 的 目标 是 求 
北海 道 庞 的 生存 数量 关 
于 时 间 的 函数 。 





那么 ， 将 北海 道 庞 的 生存 
数量 PE 随时 间 t 的 变化 做 成 
曲线 ， 则 变 成 这 样 。 


那个 …… 横 坐标 是 自 
变量 ， 纵 坐标 是 因 变 
Een 





因此 时 间 t 是 自 变量 ， 
ae ile it AR HEP 
是 因 变 量 。 


所 以 说 ， 


从 现在 开始 讲 的 
我 们 现在 首先 需要 确认 内 容 都 是 非常 重 
的 是 这 些 变 量 在 微分 方 ) 
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注意 看 这 个 微分 方程 ! 


要 目 己 想 着 做 | 然后 分 别 对 两 边 
才 行 啊 ! 进行 积分 ! 


你 小 子 还 是 钻 


钻研 什么 啊 ? | 


看 看 这 个 式 子 ， 是 不 是 
左边 集中 着 变量 已 都 到 
了 左边 ， 右 边 集中 着 变 
量 t 者 到 了 右边 ? 





西边 同时 除 以 P， 就 变 成 这 个 CN 
pp oes f 










两 个 变量 分 别 被 各 自 的 
积分 分 离 。 


这 样 的 微分 方程 叫做 可 
分 离 变量 微分 方程。 









dP 
总 结 计算 过 程 如 下 。 al aa 


eran 


二 一 一 = 儿 





BUREN, 
fxsx= tnla] 


JS zipi +c, yJat=ut+C, 


\e x 扣 职 分 常委 综合 为 一 | 
个 ， 就 变 成 放样 。 
In| P| =t+C 
由 于 PP 表示 的 是 北海 道谢 的 生存 数 
量 ， 因 此 要 注意 通常 P 是 大 于 0 的 。 


b 
微分 方程 的 解 一 




















































































































这 样 就 把 微分 方程 解 
出 来 了 。 








TVA ARR Hy BA He 


生存 数量 随时 间 的 变化 。 
Pt) =e 











那么 ， 作 图 看 一 下 ， 
会 是 什么 样 的 情况 
R? 


但 是 参数 // 和 积分 党 
数 C 都 还 没 确 定 下 来 
Pree 





无 法 具体 地 描述 这 个 
Fh, 


所 以 例 不 如 做 无 限 多 
ABA, / 


如 果 不 对 现实 世界 的 情况 进 
行 调查 ， 就 不 知道 参数 应 该 
取 什么 值 ， 积 分 常数 也 是 一 














该 常数 表示 t/0 时 北海 道 旗 
的 生存 数量 。 




















为 调查 积分 常数 ， 令 t/0。 
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像 忆 一 样 ，t=0 时 刻 的 值 
叫做 初始 条 件 。 


虽然 ， 未 必 能 够 使 用 t=0 时 刻 
的 值 。 


A 

7 =Po 

xerasiniose, 由 一 开 A 区 Pe- Fre 
ERA N PRR ) 


ATTE ELE] 1=0 时 ，P 不 随时 间 变 化 ， 
保持 二 0 时 刻 时 的 值 中 不 
变 。 


因此 我 们 假设 初 
始 条 件 及 和 比例 BEF MISH RECREDH AR 


P 











LELET LER 


e iii 指数 部 分 是 0， 因 此 不 随 


时 间 发 生变 化 。 





Jk>0 的 情况 下 ，P 随 
时 间 急 剧 增加 。 





当 / 值 变 大 时 ， 以 更 快 的 
ET Uae 





FRRONDOU SHAE CR EDES 





当 / 值 变 小 时 ， 已 以 更 
快 的 速度 减 小 。 













































































这 就 相当 于 北海 道谢 的 生 
在 数 量 连续 增加 的 情况 。 eZ 

















RB 
































作 图 能 够 更 容易 地 掌握 全 局 
变化 的 情况 。 


作 图 要 抓 住 特点 。 

















现在 的 情况 是 这 


TOG E i 
重 的 增加 








P(t)= Pye 











那么 ， 从 现实 中 北海 道 
RURE, Athik 
A Pete BRE, 


然后 推算 并 预测 一 下 今后 
北海 首席 的 数量 。 





这 是 北海 道 庞 的 估计 


值 数 据 。 


La 

















从 这 个 结果 可 以 
求 出 ， 


10.098, 
































RREBBHE E HLS BS PAE ISR 
MEBER AGE 

id 

pti 

© 
250 
200 o e 
o 

180. z 

100- be oo o o © 

=| 0°0 

o4 aR 

Assa 1956 1956 1998 2000 2002 2004 2006 2008 
根据 《关于 北海 坊 鹿 个体 指数 的 推 稼 》 的 记录 
出 作 。 注 意 这 里 的 生存 数量 是 首 定 值 ， 因 此 有 
一 定 的 误差 


另外 ， 可 以 推断 出 
2000 年 时 可 以 达到 
6 万 头 左右 ， 


因此 ， 初 始 条 件 的 t 时 
刘 定 为 2000 年 。 
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如 果 以 这 样 的 发 展 过 
度 持续 增加 ， 可 以 预 
测 出 大 概 7 年 左右 就 























[海江 西部 地 区 的 预测 
dade See Lea sin 2050 年 变 成 798 万 ，2100 
年 就 变 成 11 亿 。 
AE A MERIR 
aooo aero 
ae 2020 年 就 变 成 ‘diac i" 
万 
ooooo ao 
verano xo 
seo0000 有 a a 
。 hig 
Wo wo Ro wo Wo wo gi wo mo mo mo HO 





f 不 要 小 看 指数 
BR, 











如 果 这 个 模型 能 够 解 
释 这 个 现象 的 话 ， 





































































北海 阐 西 部 地 区 北海 并 鹿 生存 数量 的 预测 















RHZERM Maga 
pey onto 


25000000| 


20000000| 


后 面 发 生 的 现象 和 条 件 也 
会 变化 ， 因 此 也 有 必要 再 
次 探讨 模型 和 解 的 解释 。 


假设 就 算 没有 干预 
手段 ， 最 后 北海 道 
庭 的 数量 也 会 偏离 
指数 函数 发 展 。 


2000 年 的 生存 数量 为 Pu=6 万 头 ，jF=o.o8 
















当然 不 会 按 这 个 模型 









人 们 应 该 有 效 地 利用 
预测 干预 等 手段 使 / 













过 一 段 时 间 会 稳定 
在 固定 的 数量 上 。 


看 


FA 
REEN y 
的 王国 中 7 a 
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不 管 怎样 ， 需要 验证 从 模型 预测 到 的 





结果 能 否 很 好 地 解释 现实 
也 不 可 能 预测 中 的 现象 
一 次 就 解决 问 
题 。 一 让 、 
J 遇 到 不 适用 的 状况 时 ， 
























































































有 必要 进行 再 探讨 。 





















































这 部 分 比较 难 ， 再 看 
看 别 的 例子 吧 ! 


| 5. 马 尔 萨 斯 法 则 | 


© 


seo ec 


SO 89 


把 现实 世界 的 现象 模型 化 并 带 入 数学 世界 时 ， 人 们 采取 了 算式 这 样 一 种 抽象 的 表 
现 方式 。 令 北海 道 庵 的 生存 数量 为 P， 时 间 为 f， 则 通过 微分 方程 


aP 
up 
a 


表现 了 北海 道 鹿 生 存 数 量 的 增加 率 dP/di 与 北海 道 鹿 生 存 数量 P 成 正比 的 模型 。 


是 啊 ! 


但 是 ， 在 抽象 化 的 数学 世界 里 ， 用 P 表 示 北 海道 鹿 的 生存 数量 ， 完 全 不 能 让 人 产生 
兴趣 呀 。 


总 感觉 既 没 有 乐趣 也 不 讨 人 喜欢 …… 

根本 不 是 那样 的 ! 总 之 ,若是 完全 适合 的 现象 ， 用 同一 个 微分 方程 表现 的 模型 也 
可 以 解释 P 到 底 表示 的 是 什么 。 现 在 这 种 情况 下 ， 若 dP/dt 与 P 成 正比 ，P 表 示 什 么 
都 可 以 。 


“世界 人 口 ” 和 “热门 商品 的 普及 率 ” 是 这 样 吗 ? 


“细菌 的 增殖 过 程 ”也 是 。 任 何 一 个 数 的 增加 率 与 这 个 数 本 身 成 正比 ， 能 够 用 这 
个 说 法 解释 的 现象 有 很 多 。 


世界 人 口 持续 增加 ， 这 个 现象 被 人 们 所 熟知 。 
世界 人 口 大 爆发 等 …… 
实际 上 ， 成 为 问题 的 并 不 是 人 口 的 持续 增加 ， 而 是 人 口 的 指数 函数 式 增加 。 极 端 


点 说 的 话 ， 就 算 世 界 人 口 只 是 指数 函数 式 增加 ， 也 并 不 会 成 为 问题 。 成 为 问题 的 
是 人 类 生存 所 需要 的 资源 并 没有 指数 函数 式 地 增加 。 
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PIO) z 
6 
5 
4 
3 we 
2 
1 


1940 1950 1960 1970 1980 1990 2000 2010 
(公元 ) 


世界 人 口 正在 旺 指数 函数 式 增加 
令 世 界 的 人 口 ! 


© ARRERA. 


(3) 是 啊 。 比 如 说 ， 人 不 吃 东西 是 无 法 生存 的 ， 但 是 食物 的 生产 并 没有 呈 指 数 函数 式 
增加 。 若 要 使 食物 生产 像 人 口 增加 一 样 指数 函数 式 地 增加 ， 则 需要 增加 好 几 们 的 
时 间 。 一 方面 ， 通 过 增加 旱田 和 牧场 的 面积 等 食物 增产 方法 ， 能 够 提高 生产 效 
率 ， 但 是 成 倍 地 增加 旱田 或 者 牧场 的 面积 是 很 难 办 到 的 。 


© ween 


© 即使 拼命 努力 ， 在 某 个 时 期 内 也 只 能 增加 到 一 定 的 程度 ?。 成 倍 地 提高 生产 效率 是 
更 难 的 事情 。 也 就 是 说 ， 人 口 的 增长 方式 和 食物 的 增加 方式 有 本 质 的 不 同 ， 减 少 
的 不 只 是 每 个 人 的 食物 ， 而 是 随时 间 的 推 稀 ， 会 急剧 地 减少 3。 


©) 这样 下 去 ， 横 坚 都 要 陷入 饥 俄 状态 4。 
© ”最 和 提出 人 口 增长 方式 与 食物 增加 方式 不 同 之 处 的 人 就 是 马尔 萨 斯。 
©) BREE 


1 数据 出 处 : 联合 国 ， 世 界 人 口 展望 ，2008 修 订 版 。 
ee t aaa 但 是 能 够 利用 的 土地 实在 是 出 手 意料 地 少 。 另 外 ， 即 使 有 土地 ， 也 会 有 水 、 日 照 、 气 温 
面 的 问题 。 
3 也 有 等 比 级 数 式 变化 和 等 差 比 数 式 变化 。 成 售 增加 的 是 等 比 级 数 式 的 变化 ,增加 恒定 量 的 是 等 差 级 数 式 的 变化 - 
4 不 是 哪 一 个 区 域 ， 而 是 根据 地 域 的 不 同 或 许 已 经 预先 设 定好 了 。 看 看 现在 世界 的 状况 就 可 以 了 解 到 ， 不 是 全 世界 都 同时 进 
人 饥饿 状态 ， 而 是 从 贫穷 的 区 域 开始 进入 饥饿 状态 - 


现实 世界 | 数学 世界 





O 世界 人 口 增加 的 现象 


© 马尔 萨 斯 法 则 的 内 容 为 : 人 口 的 增加 率 dPyd 与 人 口 P 成 正比 。 这 个 模型 是 ， 令 时 
间 为 f 人口 为 P)， 那 么 ， 这 个 模型 是 人 口 的 增加 率 dP/d 与 人 口 PU 成 正比 。 比 例 
常数 为 上 ， 则 微分 方程 为 


SP ap 二 表示 人 口 增加 的 微分 方程 


参数 几 叫 做 增殖 率 ， 也 叫 马尔 萨 斯 径 数 '。 


© 像 细菌 这 类 微生物 也 是 这 样 ， 在 单位 时 间 里 个 体 的 增加 率 dP/dr 与 个 体 数 P 成 正 
比 。 写 成 微分 方程 的 话 ， 仍 然 是 


z. HP RIRA A EE 


E 细 疯 的 增加 速度 好 像 非常 地 快 。 


© 当然 ， 与 人 口 相 比 ， 细 菌 的 数量 可 以 说 瞬间 增加 ， 因 此 细菌 的 增殖 率 /上 与 人 口 相 
比 是 非常 大 的 数值 。 


5 注意 径 数 是 参数 ， 而 不 是 系数 。 
6 ne ， 在 条 件 适合 的 情况 下 20 ~ 30 分 钟 的 时 间 内 就 会 成 倍 地 增加 。 这 个 时 间 叫 做 世代 时 
间 。 
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完全 不 同 的 现象 篆 后 
的 法 则 确 是 相同 的 ! 


虽然 用 同一 个 模型 表示 ， 也 
不 会 有 不 协调 的 感觉 。 





表示 磁场 的 磁力 线 和 表示 液体 
流动 的 流 线 也 可 以 用 相同 的 微 
分 方程 表达 ， 








磁力 线 的 例子 














最 重要 的 内 容 ， 因 此 自己 反复 
练习 求解 才 可 以 哟 ! 





€ 
* 
= 
x 
7 
kad 
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OR 象 


北海 道 有 北海 道 鹿 等 众多 野生 动物 ， 但 是 众所周知 ， 人 类 生存 的 历史 也 较 长 。 各 
地 都 有 遗址 被 发 现 ， 其 中 甚至 有 一 万 年 前 旧 石 器 时 代 的 遗物 。 但 是 ， 如 何 才能 知道 这 些 
遗迹 是 一 万 年 以 前 甚至 更 久 的 东西 呢 ? 考古 学 中 ， 虽 然 可 以 通过 积累 地 层 的 顺序 ( 地 层 
序 ) 、 树 木 的 年 轮 等 多 种 证 据 推断 年 代 ， 但 在 这 些 方法 中 最 强 有 力 的 方法 是 利用 放射 性 
物质 的 衰变 速度 来 测定 年 代 '。 

提起 放射 性 物质 ， 总 让 人 感觉 存在 着 一 股 危 险 的 气息 ， 但 是 实际 上 自然 界 也 有 恒定 
比例 的 放射 性 物质 存在 ,平常 在 我 们 的 身边 也 存在 ?。 年 代 测定 中 ， 经 常 使 用 的 物质 是 碳 
元 素 ， 它 是 组 成 包括 我 们 人 类 在 内 的 所 有 生物 身体 的 基本 元 素 之 一 。 

形成 我 们 身体 的 碳 元 素 追 根 究 底 是 植物 通过 光合 作用 从 空气 中 的 二 氧化 碳 吸取 的 。 
植物 在 吸取 碳 元 素 的 时 候 ， 如 果 释 放出 与 时 钟 相 关联 的 东西 ， 就 可 以 得 知 这 个 植物 生存 
的 年 代 。 另 外 ,动物 的 生存 要 靠 食用 吸取 碳 元 素 的 植物 ， 也 就 是 说 ， 可 以 大 概 推测 食用 
这 些 植物 的 动物 (还 有 食用 这 些 动物 的 动物 ) 生存 的 年 代 ， 但 是 在 这 些 碳 元 素 中 ， 能 够 
那么 巧合 地 发 现 类 似 时 钟 的 东西 吗 ? 答案 是 肯定 的 。 

作为 碳 元 素 同位 素 ;的 C+，"#C 是 宇宙 中 发 射 的 放射 线 与 上 层 空气 发 生 磁 挤 所 产生 的 
中 性 粒子 ， 与 空气 中 的 氮 元 素 'N 反 应 的 产物 。 产 生 的 碳 元 素 *C 与 氧气 ( 0, ) 发 生化 合 反 
应 变 成 二 氧化 碳 ( CO, ) ， 在 大 气 中 扩散 。 

顺便 说 一 句 ， 这 个 “%C 不 能 长 时 间 稳定 地 存在 。 最 终 会 发 射 放射 线 ， 变 成 其 他 的 元 
素 ;。 不 仅 限于 "C， 发 射 放射 线 变 成 其 他 元 素 的 现象 叫做 核 衰变 ， 发 生 核 衰变 的 同位 素 叫 
做 放射 性 同位 素 。 放 射 性 同位 素 按照 恒定 的 概率 进行 训 变 。 由 于 发 生 核 误 变 变 成 其 他 元 
素 ， 所 以 原来 的 放射 性 同位 素 的 数量 ( 若 不 重新 供给 的 话 ) 会 减少 。 


1 除 此 之 外 还 有 热 冷光 法 、 电 子 回旋 共鸣 法 等 ， 但 是 考虑 方法 基本 相同 - 

2 当然 也 有 人 工 的 放射 性 物质 。 另 外 ， 即 使 是 原本 在 自然 界 中 存在 的 放射 性 物质 ， 含 量 大 时 也 会 有 危险 。 不 仅 限于 放射 性 物 
质 ， 即 使 是 存在 于 自然 界 中 的 化 学 物质 也 不 能 说 是 安全 的 。 

3 相同 元 素 的 质量 数 不 同 的 原子 ， 也 叫 放射 性 同位 素 。 原 子 核 中 的 质子 数 相同 ， 中 子 的 数量 不 同 。 

4 质量 数 为 14 的 碳 元 素 。 常 见 的 碳 元 素 是 3C， 质 量 数 为 12。 顺 便 说 一 下 ， 碳 元 素 的 原子 序数 为 6。 

5 发射 B 射线 (高 能 量 电子 ) ERN. 


原 有 的 放射 性 同位 素 的 原子 经 过 衰变 之 后 变 成 原来 一 半 的 时 间 间 隔 叫 做 半衰期 ， 
根据 放射 性 同位 素 的 不 同 而 数值 不 同 。"C 的 半衰期 是 5730 年 。 即 使 4C 与 ,发 生化 学 反应 
变 成 CO,， 被 植物 和 动物 摄取 ( 这 也 是 化 学 反应 ) ， 也 会 发 生 核 训 变 ， 若 不 提供 新 的 碳 
元 素 ，MC 经 过 5730 年 之 后 就 变 成 原来 的 一 半 。 

但 是 ， 就 像 刚才 所 说 ， 地 球 上 的 %C 是 从 上 层 空 气 经 常 性 地 供给 ， 而 经 过 误 变 数量 减 
少 的 话 ， 供 应 的 数量 就 会 相应 地 增加 以 保持 一 定 的 比例 ， 因 此 大 气 中 ( 以 CO, 的 形式 ) 
存在 的 *C 的 量 是 大 致 不 变 的 ?"。 也 就 是 说 ， 大 气 中 的 *C 和 "3C 的 比例 是 恒定 的 *。 

由 于 含 %C 的 CO 和 含 C 的 CO, 的 化 学 性 质 完全 相同 *， 因 此 无 区 别 地 被 植物 吸收 进行 
光合 作用 。 正 在 进行 光合 作用 的 植物 细胞 所 含 的 2C 和 "*C 之 比 与 大 气 中 的 比例 应 该 是 相同 
的 。 也 就 是 说 ， 植 物 细胞 在 进行 光合 作用 的 时 候 ， 其 细胞 内 的 2C 与 AC 的 比例 也 是 恒定 
的 。 

但 是 ， 当 植物 细胞 不 进行 光合 作用 的 时 候 ， 就 没有 大 气 给 植物 提供 8AC， 因 此 '2C 对 
#C 的 比例 会 减 小 。 也 就 是 说 ， 当 植物 不 能 吸收 碳 元 素 的 时 候 ， 就 会 放出 类 似 时 钟 的 东 
西 。 不 能 进行 光合 作用 的 理由 可 以 考虑 为 枯萎 、 变 成 木质 %*、 被 动物 吃 掉 等 。 不 管 是 动物 
还 是 植物 ， 通 过 调查 生物 遗体 内 所 含 ?C 与 *C 之 比 就 可 以 得 知 光 合作 用 是 何 时 停止 的 。 


6 有 些 放射 性 同位 素 的 半 误 期 仅仅 是 一 瞬间 ， 而 占 天 然 钠 元 素 大 部 分 的 钠 238(U) 的 半 误 期 有 45 亿 年 之 长 。 顺 便 说 一 下 ,用 
于 核电 和 核武 器 的 铀 235(3U) 的 半衰期 是 7 亿 年 。 
be 如 果 速 度 相 同 ， 就 始终 处 于 同一 高 度 ， 与 这 件 事情 是 相同 的 道理 。 但 是 好 孩子 是 不 会 那样 模仿 


a emu NEE, 为 此 ， 年 代 测定 的 正确 性 受到 质疑 ， 关 于 推测 年 代 的 讨论 也 是 此 起 彼 伏 。 


9 物理 性 质 不 同 。 比 如 说 质量 是 中 子 的 两 个 大 小 ， 因 此 扩散 速度 要 慢 一 些 。 
10 越 车 近 树木 年 轮 的 内 侧 "C 的 比例 就 越 小 。 
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= t 
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直到 光合 作用 结束 为 止 是 恒定 的 ， 但 是 光合 作用 结束 之 后 每 过 5730 年 就 会 减少 为 原 
来 的 一 半 。 相 反 ， 若 得 知 "C 与 “C 之 比 ， 就 能 得 知 光合 作用 是 何 时 停止 的 。 


仿生 物体 内 所 含 *C 与 *C 之 比 





虽然 刚才 说 到 ， 放 射 性 同位 素 按 照 恒定 的 概率 进行 衰变 ， 但 是 如 果 只 观察 某 个 放射 
性 同位 素 的 一 个 原子 ， 说 实话 现在 还 不 能 得 知 这 个 原子 何 时 开始 进行 核 训 变 。 虽 然 现在 
不 行 ， 也 许 在 更 遥远 的 未 来 有 可 能 得 知 。 某 个 原子 进行 衰变 的 时 机 是 随机 的 。 不 管 是 与 
uC 种 类 相同 的 放射 性 同位 素 中 的 哪个 原子 ， 若 继续 观察 ， 都 是 以 相同 的 概率 进行 衰变 ， 
这 一 点 是 可 以 确定 的 。 即 使 仅仅 观察 一 个 原子 ， 也 不 能 预测 该 原子 从 何 时 开始 衰变 的 
但 是 当 观 察 很 多 个 原子 的 情况 下 ， 可 以 预测 到 随 着 时 间 的 推移 会 有 恒定 比例 的 原子 开始 
衰变 。 也 就 是 说 ， 原 子 数 的 变化 率 ( 衰变 速度 ) 与 原子 数 成 正比 。 认 ? 好 像 是 在 哪儿 见 
过 的 词汇 啊 。 对 了 ， 是 漫画 页 中 对 北海 道 鹿 的 生存 数量 P(D) 的 增加 率 dP(D/dx 与 北海 道 鹿 
的 生存 数量 Pi) 成 正比 的 模型 进行 的 假设 。 这 里 也 要 相同 地 把 核 训 变 模型 化 之 后 带 人 数学 
世界 。 
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放射 性 同位 素 衰变 现象 的 说 明 。 


O 核 衰变 的 说 明 


ge 型 

那么 ， 试 着 做 一 个 说 明 放射 性 物质 误 变 的 模型 。 模 型 化 的 时 候 需 要 留意 的 是 即使 不 
知道 这 种 现象 为 什么 发 生 也 没有 关系 "。 当 然 ， 想 要 知道 为 什么 不 能 预测 某 个 原子 从 何 时 
开始 衰变 的 机 理 的 这 种 心情 是 可 以 理解 的 。 知 道 的 感觉 肯定 要 比 不 知道 时 更 加 顺畅 。 但 
是 ， 有 很 多 现象 虽然 不 知道 为 什么 发 生 ， 但 能 够 知道 是 如 何 发 生 的 。 暂 时 把 机 理 的 阐明 
搁置 起 来 ， 仅 仅 试 着 做 一 下 现象 的 说 明 。 从 现在 开始 应 该 能 够 明白 很 多 的 事情 了 。 

进行 现象 的 说 明 时 ， 需 要 理解 一 下 这 个 现象 是 如 何 发 生 的， 模型 化 之 后 把 它 巧妙 地 
代入 数学 世界 就 足够 了 。 虽 然 不 知道 核 衰变 现象 通过 什么 样 的 机 理发 生 ， 但 得 知 原子 数 
的 变化 率 与 原子 数 成 正比 之 后 ， 就 能 够 把 它 带 入 数学 世界 了 。 

令 时 间 为 !、 放 射 性 物质 的 原子 数 为 Nr)， 假 设 原子 数 的 变化 率 dN(1)/ dy 与 原子 数 MD) 
成 正比 的 模型 。 比 例 常 数 为 A ， 则 


aN 


go 描述 原子 数 随 核 衰变 变化 的 微分 方程 


当 原 子 数 N 随 时 间 减 小 ， 比 例 常数 入 变 成 正 值 时 ， 要 加 上 负 号 ?。 与 考虑 北海 道 鹿 生 
存 数量 时 的 情况 相同 ， 仅 仅 从 这 个 模型 不 能 决定 参数 和 3。 





| 
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由 核 误 变 产生 的 原子 数 的 变化 率 与 原子 数 成 正比 ， 因 此 得 到 原子 数 的 变化 率 duMUIJd' 与 原子 数 NU1) 成 正比 的 
微分 方程 的 模型 


O 说 明 放射 性 同位 素 豪 变现 象 的 模型 


11 不 仅 限于 核 衰变 ， 也 适用 于 一 般 情况 。 建 立 北海 道 麻生 存 数量 模型 时 ， 也 说 过 同样 的 话 。 
12 虽然 参数 本 身 也 可 以 是 负 的 ， 但 是 这 里 要 给 人 一 种 减少 的 感觉 - 
13 就 像 后 面 提 到 的 那样 ， 放 射 性 同位 素 根据 种 类 不 同 误 变 的 速度 也 会 不 同 ， 与 此 相对 应 的 入 值 也 会 不 同 。 
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那么 ， 求 一 下 这 个 微分 方程 的 解 。 把 77 页 求 出 的 式 子 AE = Pst (3.1) 相 比 之 后 发 
现 ， 形 式 几乎 完全 相同 “。 若 微分 方程 的 形式 相同 ， 则 解 的 形式 也 会 相同 。 也 就 是 说 ， 已 
经 知道 其 解 了 5， 在 这 里 就 当 练习 一 下 求解。 

我 们 的 目标 是 求 放 射 性 物质 的 原子 数 关 于 时 间 的 函数 Ntt)。 每 次 遇见 微分 方程 ， 都 
需要 确认 一 下 微分 方程 的 这 些 变量 。 


自 变量 


两 边 再 进行 积分 ， 则 


1 
fav- -afa 
左边 的 变量 N 和 右边 的 变量 t 分 别 被 各 自 的 积分 分 离 。 


1 
nes 
求 各 自 的 积分 得 到 
ax=ajwlrc 
-afa =-M+0, 
综合 一 下 积分 常数 得 到 


In|N|=-\t+C 


14 意 恩 是 ， 仅 仅 做 变量 替换 。 
15 相当 于 85 页 图 形 中 的 w<0 的 情况 。 


100 


放射 性 物质 的 原子 数 关于 时 间 的 函数 Ntt) 为 “ 
NE) =e} EAE (3.2) 
这 样 就 求 出 了 微分 方程 的 解 。 








现实 世界 | 数学 世界 N 4 
求解 微分 方程 之 后 ， 得 知 放射 性 物质 的 原子 数 随时 间 的 变化 可 以 用 指数 函数 来 表达 。 
O 说 明 放 射 性 同位 素 认 变现 象 的 模型 的 解 

af 释 

其 次 ， 就 是 决定 积分 常数 ，1=0 时 刻 的 放射 性 物质 的 原子 数 设 为 M0)=N,， 那 么 以 此 
为 初始 条 件 ， 根 据 式 (3.2 ) 则 

MO)=er 

因此 化 为 

Ne 

表示 放射 性 物质 原子 数 的 函数 能 够 写成 

ND)=Ne™ (3.3) 

由 于 已 经 决定 了 积分 常数 ， 再 假设 参数 入 ， 作 为 放射 性 物质 的 原子 数 N 随 时 间 ! 变 化 
的 曲线 。 由 于 现在 的 情况 下 入 >0， 原 子 数 必然 会 减少 。 


16 N 为 放射 性 物质 的 原子 数 ， 因 此 注意 N 恒 大 于 0。 
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参数 和 越 大 ， 衰变 越 快 。 这 里 作出 了 与 三 个 相对 应 的 三 条 曲线 。 
人 表示 放射 性 物质 的 原子 数 N 随 时 间 ! 变 化 的 曲线 


参数 \ 越 大 ， 原 子 的 数量 减少 越 快 (图 中 的 曲线 画 出 了 与 三 个 \ 相 对 应 的 三 条 曲 
线 ) 。 根 据 放射 性 同位 素 的 不 同 ， 误 变 的 速度 也 会 不 同 ， 误 变 快 的 放射 性 同位 素 的 入 值 
BRK, 衰变 缓慢 的 放射 性 同位 素 的 值 较 小 。 为 了 表示 衰变 的 程度 ， 叫做 衰变 常数 。 
研究 一 下 实际 的 放射 性 同位 素 并 得 到 其 衰变 常数 ， 则 能 够 得 知 何 时 剩余 多 少 原子 ， 也 能 
够 说 明 未 来 和 过 去 的 现象 。 
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得 到 了 衰变 常数 ， 便 能 够 得 知 何 时 剩余 多 少 原子 ， 也 能 够 预测 未 来 的 情况 。 
O 说 明 放射 性 同位 素 衰变 现象 的 模型 的 解 的 解释 


我 们 也 试 着 研究 一 下 半衰期 与 衰变 常数 的 关系 。 把 半衰期 写成 ha，t=0 时 刻 的 放射 
性 物质 的 原子 数 为 NW，:=#o 时 刻 ， 变 成 YN2， 那 么 根据 式 (3.3 ) 
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得 到 


Ni yga 
2 o 


SRT EAA, TURE EER, ON 


| 
+ -了 人 | 0.69 


本 Caan © 

也 就 是 说 ， 半 误 期 与 衰变 常数 的 倒数 成 正比 "。 

那么 ， 我 们 能 够 测定 的 是 放射 性 强度 *。 放 射 性 强度 与 单位 时 间 内 衰变 的 原子 数 
(原子 的 变化 率 ) 成 正比 ， 因 此 也 就 是 与 原子 数 成 正比 5*。 需 要 说 明 的 是 ， 豪 变 常数 不 能 
直接 进行 测定 得 到 ， 而 是 根据 放射 性 强度 随时 间 减 少 的 情况 求 得 。 

哎呀 ， 我 们 终于 理解 了 如 何 利用 放射 性 物质 测定 年 代 。 放 射 性 物质 发 射 放射 线 的 同 
时 误 变 成 其 他 的 元 素 ， 我 们 了 解 到 在 这 个 过 程 中 放射 能 的 变化 如 式 (3.3) 所 示 ， 对 于 已 
知 半衰期 ( 也 就 是 衰变 常数 ) 的 放射 性 同位 素 ， 若 已 知 现在 的 数 N 与 其 对 应 的 过 去 某 个 
时 刻 的 数 N, 之 比 ( 放射 性 强度 之 比 ) NN,， 就 可 以 得 知 其 时 间 间 隔 。 


—s=e" 


N, 


A 
_ ml 

入 

这 真是 太 棒 了 。 也 就 是 说 ， 只 要 得 到 从 遗址 出 土 的 遗物 的 MN,， 就 能 知道 这 些 东 西 
是 什么 时 候 制作 的 了 2。 

从 遗址 出 土 的 文物 ， 比 如 出 土 的 是 木 制品 ， 就 可 以 得 知 这 些 树 何 时 停止 光合 作用 
的 。 可 是 并 不 能 得 知 树木 被 切割 的 年 代 和 被 加 工 的 年 代 ， 这 一 点 虽然 较为 玄妙 ， 但 是 那 
个 时 代 的 人 不 去 特意 挖掘 并 使 用 很 久 以 前 的 木材 ， 就 可 以 推测 出 遗迹 的 年 代 。 如 果 测定 
一 下 那 一 带 的 碎 木片 什么 的 ， 也 许 是 年 代 更 加 久远 的 东西 2。 





17 衰 变速 度 快 的 话 ， 变 成 原来 一 半 的 时 间 就 会 变 短 。 因 此 这 是 理所当然 的 。 
18 测 定时 ， 提 起 放射 能 ， 经 常用 到 的 是 盖 革 -弥勒 计数 管 和 办 光 检 出 器 。 为 了 能 够 简单 地 得 知 被 爆 线 量 而 使 用 胶片 式 射线 测 
定 器 


19 就 是 说 ， 如 果 放射 性 同位 素 的 原子 数量 较 多 的 话 ， 放 射 性 也 较 强 。 

20 正 规 地 说 ， 是 什么 时 候 停止 光合 作用 的 。 

21 这 种 根据 放射 性 同位 素 测定 年 代 的 方法 ， 年 代 越久 远 测定 难度 就 越 大 。MC 的 半 误 期 为 5730 年 ， 大 概 6000 年 左右 变 成 原来 
的 一 半 ,，12000 年 就 减少 为 原来 的 四 分 之 一 ， 过 了 60000 年 之 后 就 变 成 原来 的 11000。 量 变 的 少 , 就 很 难 进行 正确 的 测定 。 
能 够 用 "C 进 行 测定 的 也 只 有 数 万 年 的 限度 。 超 过 此 范围 ， 就 不 能 够 很 准确 地 测定 年 代 更 久远 的 东西 了 。 若是 岩石 ， 有 方 
法 可 以 测定 年 代 更 加 久远 的 东西 。 
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E 各 种 各 入 的 现 旬 与 一个 表达 式 ] 








到 此 为 止 ， 我们 得 知 说 明 马尔 萨 斯 法 则 与 核 衰变 这 两 个 完全 不 同 现象 的 微分 方程 在 
数学 上 可 以 用 一 个 形式 表达 。 到 现在 为 止 解决 的 可 分 离 变 量 微分 方程 在 微分 方程 中 是 最 
基本 的 。 方 程 中 所 含 的 两 个 变量 被 各 自 的 积分 分 离 ， 能 够 进行 积分 ， 就 可 以 解 出 这 个 微 
分 方程 。 解 法 虽然 简单 ， 但 是 各 种 各 样 的 现象 都 能 用 这 个 形式 的 微分 方程 来 表达 。 

另外 ， 我 们 得 知 说 明 物体 冷却 时 的 温度 随时 间 的 变化 ( 牛顿 冷却 法 则 ， 参 考 附录 
1) ,计算 火箭 的 到 达 速度 ( 齐 奥 尔 科 夫 斯 基 公 式 ， 参考 附 录 2 ) 对 应 刺激 的 心理 量 的 变化 
( 费 希 纳 的 法 则 ， 参 考 附录 3 ) 等 各 种 各 样 的 现象 都 能 用 可 分 离 变量 微分 方程 表达 。 数 学 
世界 中 形式 相同 的 模型 ， 能 够 说 明 现实 世界 中 各 种 各 样 的 现象 ， 是 不 是 有 点 神秘 的 感觉 
w? 
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L .物流 模型 | 


那么 ， 在 91 页 研究 马尔 萨 斯 法 则 时 ， 使 用 了 世界 人 口 持续 增加 这 个 话题 。 但 是 ， 最 
近 实际 人 口 的 增加 率 有 所 下 降 。 


POIO) i 


1940 1960 1980 2000 2020 2040 2060 
{公元 ) 


© 世界 的 人 口 (预测) ? 
一 般 来 说 ， 如 果 人 口 一 味 地 呈 指 数 函数 式 增加 ， 不 久之 后 世界 人 口 的 质量 就 会 超过 
地 球 本 身 的 质量 。 不 管 怎 么 样 ， 那 样 的 事情 总 不 会 发 生 的 。 马 尔 萨 斯 法 则 虽然 局 部 地 说 
明了 一 些 现象 ， 但 是 逐渐 与 现实 发 生 了 偏离 。 














Ea 现实 世界 | 数学 世界 


虽然 ， 之 前 说 明 的 人 口 增加 现象 中 提 到 增加 率 最 近 有 所 下 降 ， 但 是 解 油 分 方程 得 到 的 结果 是 增加 率 在 持续 上 升 。 由 于 对 这 样 
的 结果 不 满意 ， 因 此 再 一 次 返回 到 数学 世界 对 模型 进行 修正 。 


© 人 口 增加 模型 的 修正 
22 数据 出 处 : 联合 国 ， 世 界 人 口 展望 ，2008 年 版 。 
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不 仅 限于 人 口 。 虽 然 刚 开始 的 时 候 ， 皮 式 培养 吗 中 培养 的 细菌 也 是 按 指数 函数 式 地 
进行 增殖 ， 但 是 数量 增加 到 不 能 从 培养 池 里 吸收 营养 的 时 候 就 停止 增殖 ， 增 加 率 的 变化 
就 会 迟缓 >。 在 稍微 大 一 点 的 范围 使 用 这 个 模型 时 ， 就 不 得 不 对 这 个 模型 进行 修正 。 

想 说 明 某 个 现象 的 时 候 ， 将 现实 世界 的 现象 带 和 数学 世界 ， 然 后 试 着 再 用 数学 世界 
里 解 出 的 结果 解释 并 说 明 这 个 现象 。 尽 管 得 到 了 期 盼 的 结果 并 且 能 用 它 结束 这 个 作业 ， 
但 是 在 不 够 满意 的 情况 下 ， 再 一 次 回 到 数学 世界 对 模型 进行 修正 ， 反 复 求解 解释 这 样 的 
作业 。 这 样 一 来 ， 慢 慢 地 就 能 得 到 更 加 精密 的 解 。 要 精密 到 什么 程度 才 行 呢 ? 到 能 够 认 
可 的 程度 就 可 以 了 。 进 行 循环 之 后 ， 比 以 前 还 要 更 上 一 层 ， 因 此 与 目标 的 距离 就 逐渐 缩 
短 了 。 那 样 还 不 能 达到 满意 的 话 ， 再 一 次 回 到 循环 中 ， 对 剩 下 的 那 一 些 内 容 进行 进一步 
的 研究 和 探讨 。 

再 考虑 一 下 人 口 增加 的 情况 。 人 口 增加 ， 仅 仅 对 生活 环境 的 影响 也 是 很 大 的 *。 
生活 环境 恶化 的 话 ( 越 来 越 难以 居住 ) ， 增 加 率 也 会 有 下 降 趋势 。 这 里 ， 把 马尔 萨 斯 径 
数 改 写成 随 着 人 口 P 的 增加 变 小 的 形式 。 设 K 为 常数 ， 则 改写 为 


HB) 
则 微分 方程 变 成 
dP 


p? ` 
P -zl- 天) P。 -被 修正 过 的 措 述 人 口 增长 的 微分 方程 (34) 


那么 ， 这 个 微分 方程 的 解 曲 线 会 是 什么 样 呢 ? 
整理 式 (3.4 ) ， 右 边 变 成 以 K 为 分 母 的 形式 ， 则 


aP _, (K-P)P 
a M K 
虽然 有 点 难 懂 ， 但 还 是 可 分 离 变 量 微分 方程 。 
BB_ KO 
dt K 
变量 分 离 
a: -上 
je- 和 i 


其 次 ， 要 对 其 进行 积分 ， 右 边 倒是 没有 问题 ， 左 边 的 话 就 有 点 难 办 了 。 这 时 会 用 到 
一 些 常用 的 便利 小 招数 。 将 左边 的 积分 函数 改写 为 


23 增殖 好 像 可 以 分 为 好 几 个 阶段 
24 环境 问题 。 








(B+ sa) ) 
(K-PP K\P K-P 
进行 通 分 的 逆 运 算 。 分 解 部 分 分 数 。 
a({a 1 
(| Laps xe dj- 中 
那么 ， 到 此 为 止 各 项 的 积分 都 能 用 我 们 前 面 见 过 的 方法 进行 处 理 。 对 括号 里 的 两 项 
积分 ， 进 行 变量 替换 ， 令 s=K-P， 根 据 ds/dP=-1 则 
1 1dP 1 
[ph par-J 482 ana | tare tneina-P) 
(这 个 叫做 变量 替换 法 ) ， 整 个 式 子 化 为 *” 


InP-In(K- P)=pt+C 
利用 对 数 的 性 质 


P 
ne P 

写成 指数 函数 形式 
Pane 
Kp 

求解 P， 则 能 够 求 出 

Kee 

P(t)= ae 全 被 修正 的 微分 方程 的 解 (3.5) 

设 /=0 时 刻 的 人 口 为 P,， 则 

__Kee 

”1+ec 


由 上 式 得 到 





= pt +0 








把 上 式 代入 式 (3.5) 中 进行 整理 
KPe* 
K+ P(e*-1) 
得 到 的 解 虽然 有 点 复杂 ， 这 时 ， 要 领 就 是 从 极端 的 地 方 开始 考虑 。 


PG= (3.6) 


25 P、( K-P) 同 为 正 。 
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比如 说 ， 把 上 0 代入 式 (3.6) ， 就 变 成 P(0)= P,。 这 是 初始 条 件 中 交代 的 内 容 ， 理 所 
当然 成 立 。 试 求 :> 时 的 极限 ， 但 是 从 式 (3.6) 不 容易 看 清楚 ， 分 子 和 分 母 分 别 除 以 
em, 再 整理 

KP, 
PH= 一 一 一 
(K-P)e*+ P, 
HHR MHR, P =K WREN, RAP (0) = P, 开 始 上 升 ，P。=K 时 
应 该 变 成 最 高 曲线 。 作 出 其 解 曲线 。 


P 





5 
刚 开始 的 时 候 ， 以 指数 函数 ( 细 线 ) 式 地 增加 ， 接 着 增加 率 变 得 迟缓 ， 然 后 再 到 达 恒 定 的 值 达到 他 和 。 
O 被 修正 的 人 口 增加 模型 的 解 曲线 






po 


7 
[ 

L | 
| t 












现实 世界 | 数学 世界 


进行 了 两 个 循环 之 后 得 到 了 期 望 的 结果 。 
依 被 修正 后 的 人 口 增加 模型 


这 样 的 模型 被 称 作物 流 模型 。 不 仅仅 适用 于 人 口 增 加 ， 也 能 够 很 巧妙 地 说 明 工业 
制品 普及 的 情况 。 

本 章 主要 讨论 了 可 分 离 变量 微分 方程 。 这 里 大 体 地 看 一 下 

一 阶 微分 方程 


p(s)G(y) (38) 
dy 


中 ， 右 边 是 z 的 函数 F(z ) 与 y 的 函数 G (y) 的 乘积 。 像 这 种 形式 的 微分 方程 ， 两 个 变 
量 z 和 y 分 别 能 够 被 分 离 到 左边 和 右边 ， 因 此 叫做 可 分 离 变 量 微分 方程 。 
式 (3.8 ) 的 两 边 除 以 G (y), 


两 边 同时 对 x 进行 积分 ， 


| | F(z)dz+C 
Gly) dz 


左边 的 积分 用 变量 替换 法 改写 成 


设 积分 常数 为 C， 则 化 为 


| F(z)dz+C (3.9) 

左边 是 只 含 变量 y 的 微分 dy 和 变量 y 的 函数 G (y) ， 右 边 是 只 含 变量 z 的 微分 dz 和 
变量 z 的 函数 尺 ( > ) ， 因 此 得 知 左边 和 右边 可 以 进行 变量 分 离 。 变 量 分 离 之 后 ， 仅 仅 进 
行 积分 运算 就 可 以 了 。 然 而 ， 具 体 求解 的 时 候 也 会 出 现 解 不 出 来 的 情况 。 也 可 以 利用 
计算 机 进行 求解 。 将 其 变 成 式 (3.9) 的 形式 ， 一 样 能 够 解 得 出 来 x。 

引入 变量 分 离 法 也 同样 可 以 对 微分 方程 进行 求解 ， 因 此 变量 分 离 法 也 是 解 微分 方 
程 的 基本 方法 。 即 使 是 乍 一 看 不 能 看 出 是 可 分 离 变 量 微分 方程 ， 也 可 通过 巧妙 地 替换 
变量 等 招数 进行 变量 分 离 。 就 像 解 字谜 游戏 一 样 ， 虽 然 刚 开始 的 时 候 自己 做 的 话 也 许 
会 觉得 很 有 意思 ， 但 还 是 先 学 习 一 下 先 人 的 智慧 吧 。 


26 也 有 这 样 求解 的 情况 。 
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人 都 说 自己 是 非 专业 的 ) 
+, 竟然 …… A 
































































































































RABEWRA, $4 øA 

者 肯定 更 少 ， 真 是 越 来 有 什么 好 高 兴 的 ! A 

越 闲 …… an \ 
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HEET 
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四 

A 我 还 能 很 能 
于 的 | 
Kaaa Y 


Eo 
wR 
OE 
KE 
Mie 
+R 
RF 





你 们 的 愿望 不 是 已 经 


实现 了 吗 ! ? 
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Ce 


ath 
还 是 先 学 习 吧 ! P 
ai 





在 附近 的 面包 屋 的 所 有 东西 
里 ， 这 个 是 最 好 吃 的 了 ! 



























































































在 此 之 前 我 们 研究 了 递增 
递减 的 模型 。 











在 现实 中 持续 增加 持续 
减少 的 事情 不 太 常见 。 
ey 
x 
~~ 


“oy 





因此 ,今天 我 们 挑战 一 下 
逐渐 变化 然后 变 成 恒定 什 





第 4 章 多 一 阶 非 齐 次 线性 微分 方程 常数 变易 法 115 


这 是 以 能 用 变量 分 
REMMI BATE OL 
BRAM, 现象 的 说 明 吧 ! 














要 用 “从 空中 落下 ” 
的 东西 说 明 吗 ? 





这 也 是 在 落下 









你 的 意思 是 说 ， 仅 仅 
因为 大 落下 的 速度 就 








这 是 因为 速度 
ARR. 



































相反 也 有 像 雷 阵 
雨 似 的 快速 落下 
HA, 









































AAAPHKÄHh, 
而 雷阵雨 的 雨滴 较 大 的 


























重力 





那么 从 相同 高 度 落下 
的 所 有 水 滴 都 以 相同 


空气 3 对 ,水 消 在 空气 中 | 
下 落 的 时 候 会 受到 空气 的 
BA. / 





运动 的 物体 受 拖 找 空 气 所 




















产生 的 粘性 阻力 和 


运动 的 物体 从 正面 碰撞 空气 的 
时 候 产 生 的 惯性 阻力 。 


他 们 共同 作为 阻力 
HAT HR, KA 
向 为 阻碍 物体 运动 











顺便 说 一 下 ， 吉 性 阻力 
与 物体 的 速度 成 正比 
惯性 阻力 与 速度 的 平方 
REM, 





MOHELA, RRR 水 滴 刚 开始 落下 


快 ， 阻 力 就 会 越 大 吗 ? 


随 着 迷 度 越 来 越 快 ， 
阻力 也 会 越 来 越 大 ， 


以 至 于 不 能 继续 
Wik, 














这 么 说 ， 阻 力 的 大 
小 不 可 能 超过 重力 






1 acer 

F=-aLlm - BS — pv" 
M ae 
sanpa 空气 的 密度 
PEAKS PEORRER 


i emh meamh 
— — 



































j mol 
这 里 的 a 和 B 是 ? -@L nv S 7 pu 


也 就 是 说 ， 最 终 会 达到 
一 个 恒定 的 速度 。 





由 物体 形状 决定 的 常数 。 
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ee 

水 滴 受 表面 张力 的 作用 mae. OC 那么 , L 和 5 将 
SRA AH 
HERR? 






































入 
(多 体 的 大 小 )= 27 














,9 (物体 的 模 玲 面积 )= 7077 





如 果 是 球体 ， /表示 
的 是 直径 。 





那么 ， 如 果 物 体 的 形 
所 受 的 空气 阻力 就 变 
RANE, BRAT 


就 可 以 得 知 
a =3T、B =1/2, 


首先 ， 千 性 阻力 与 示 
其 相对 应 的 惯性 阻 ran 
力 之 比 。 


我 们 把 这 里 出 现 的 
Lpvih 叫 做 雷诺 数 。 





那么 ， 由 于 大 气 中 
p=l.2kg/m’, =1.8X10°Pa's, 





这 样 一 来 , L=0.1mm 
为 丫 性 阻力 还 是 惯性 
阻力 发 挥 主要 作用 的 
分 界 点 。 














SEZRRIHHR, BEB 
力 发 挥 主要 作用 ， 大 的 时 候 惯 
性 阻力 发 挥 主要 作用 。 





122 





Adie am + @ 
大 小 如 右边 所 示 。 | | 1000pn(= 1mm) 








1O pm (= 0.01mm ) 

















结构 不 太一 样 

ai BB 弹 是 介 于 雨滴 与 平衡 球 
云 履 结核 分 界 点 的 
物质 。 人 




























































我 们 试 着 把 如 云层 结核 
大 小 的 水 满 的 降落 过 程 
模型 化 。 





首先 ， 物 体 在 空气 中 
下 落 的 运动 方程 中 ， 
以 笔直 向 下 的 方向 为 
正方 向 。 
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云 活 结核 的 情况 就 是 只 考虑 车 性 阻力 。 
这 样 。 
dv 










gy mg - 6anrv 
m—=mg- 
sai PR 
m 
dv 
dt 


6rnrv 
=g- 一 -一 
m 






RMA ZA 
运动 微分 方程 


为 了 更 容易 看 性， 把 方 
程式 中 楼 度 以 外 的 量 综 


合 为 
































这 次 只 是 稍微 
含有 9， 因 此 不 能 进行 \\ C A] 提高 了 一 点 难 


RED Boe \y \ 
NAN AA 2. 






























































































































































云 疑 结核 的 返 动 微分 方程 








MERRER) | 单位 质量 物体 上 右边 表示 作用 在 物体 











上 的 力 ( 除 以 质量 的 
结果 ) 





这 个 微分 方程 的 左边 
表示 物体 过度 随时 间 
的 变化 。 
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也 就 是 说 ， 物 体 的 速度 


在 右边 力 的 作用 下 随时 、 那么 
间 变 化 。 


由 于 改变 物体 速度 的 原 
因 是 9 和 -两 项 之 和 ， 






因此 它们 决定 了 物 
体 的 速度 。 





这 两 项 之 中 ，9 表 示 重 力 加 


速度 这 个 常数 ， Ug ~ 


-表示 与 物体 OK) 的 4 
ERRENA, 
pac 
第 一 项 的 常数 表示 
物体 总 是 受 一 定 的 





























































































































从 物理 方面 考虑 2a 
— 下 。 
oC 
VACHE 7 3) 
J 
也 就 表示 物 休 以 和 要 直 
. 向 下 的 加 速度 运动 。 





































































































































































































































































































































































































SANTEE RL 以 委 直 向 下 的 加 速度 要 直 向 下 的 速度 0 
运动 的 意思 是 


右边 的 项 表示 重力 与 
阻力 的 合力 ， 


物体 下 落 的 速度 越 快 ， 


物体 垂直 向 下 的 合力 就 


peta RRL, LAMA 


保持 平衡 时 作用 在 物 
体 上 的 合力 为 0， 


阻力 大 于 重力 的 情 
讽 是 不 存在 的 1 


就 像 列 才 提 到 的 
小 水 滴 的 情况 | 

















































































































.那么 ， 试 着 考虑 一 
T, RAEN, HAE 
成 什么 样 呢 ? 





没有 重力 的 情况 下 ， 由 于 没 
ALLAH, 


刚 开始 处 于 静止 状态 的 
物体 不 会 落下 ， 会 原 地 





J, 
` A 
J 
i / 
站 
多 从 | # | ; : 
IV e So 
; 这 个 就 叫做 思想 
实验 。 


但 是 什么 也 不 发 生 
就 没有 意思 了 。 


因此 我 们 假设 物 
体 刚 开 始 时 稍微 
运动 了 一 下 。 




















AA 


设 物体 开始 运动 时 的 初速 度 为 
vw， 刚 开始 时 作用 在 物体 的 阻 
力 为 -vo。 


RH, BAHAR 
体 运动 的 方向 ， 因 此 物体 的 
RARER ED 





最 终 会 停止 的 吧 。 用 微分 方程 表示 这 个 状况 。 a) 


没有 重力 ， 表 示 第 一 项 的 8 
没有 的 情况 。 


运动 方程 只 考虑 有 粘性 阻力 的 情况 。 





dv 
-一 = 一 IJ 
oS Ce) 
fiiwr at 
dS a) 


| lgv=Invl+ C, -f dt=-1t+C, 
v 


\ AE 
REN 
In|v|=-vt + C, 
Voom 


RA EBB HH A ( ) 


= 














10H, v(0)=vo。 


假设 刚 开始 时 物体 处 于 运动 
RS, MRD 


最 后 ， 将 其 代入 刚才 得 到 总 结 一 下 得 到 如 下 所 
的 式 子 中 ， 未 知 的 常数 就 示 的 解 。 
变 成 …… 











































































































































































































KARE, RARE 
. 由 于 有 空气 得 力 ， 若 没 
Ts 有 类 似 重力 的 外 力 ， 物 
休 最 后 就 会 停止 运动 
RRRHR 
接近 于 0 
„roenn 
RRRA Jy 


i 


顺便 说 一 下 ， 设 初始 速度 bo 
为 0， 则 从 也 能 够 表达 一 直 部 保持 静止 


ut) = ne" ER 














那么 ， 在 没有 重力 的 时 
候 ， 可 利用 分 离 变 量 微 
分 方程 求解 。 

















DS paniti- 

Cece 
WEBS 

SAt+BeRAKRY 








RARER ADE | 
z 


hs 
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那么 ， 假 设 U(D)= ce "中 的 
系数 c 为 时 间 1 的 函数 ! 















































































































































































































































那么 ， 微 分 方程 的 解 U(D= ce 
可 以 写成 U(D)= ce", 








ee 
u(t) = c(t)e 

























































满足 这 个 解 的 微分 方程 就 
是 我 们 本 来 起 要 求解 的 微 
分 方程 ， 

























































































将 此 解 代入 本 来 想 要 求 的 
方程 中 ! 





第 4 章 儿 一 阶 非 齐 次 线性 微分 方程 常数 变易 法 133 


134 





所 有 计算 归纳 总 结 如 下 所 示 








elt) ou 
dt 


de 


‘+ e(t) var 


本 来 想 要 求解 的 微分 方程 
代入 假设 的 解 xD=c(D)e- 


= BE) ny op _ve 
= etre) 








>) 消去 两 边 的 vc(jje* 


2 将 e" 移 项 ， 然 后 进行 积分 


) 代入 假设 的 解 v() = ele” 


考虑 重力 和 黏 性 阻力 时 的 解 





























用 含有 微分 方程 自 变 量 的 
RERRAK, 


SA+ReRAGRY 











RERBEBA 
时 的 现 家 
i 


— 

















HL, AR PERK 
下 落 的 情况 吧 。 


假设 初始 条 件 !=0 时 wv(0)=0， 
然后 进行 求解 。 


好 像 是 这 样 吧 ? 








初始 刍 度 为 0 的 情况 下 ， 空 气 中 下 落 i 
MAH MEL, RAR 
oA 28 AE PADI WY EM MAE TRB. 





I 





这 叫做 终端 速度 。 





该 式 综合 了 粘性 
阻力 系数 和 物体 
质量 。 





利用 这 个 ， 再 还 原 为 Dy 
Rimh Akh 
质量 时 ， > 














HETER KH 


olt) = 20-4) & 
极限 得 到 的 值 。 


终 闹 速度 与 重量 (EH 
大 小 为 mg) 成 正比 ! 


那么 存在 空气 阻力 
的 条 件 下 ， 重 的 东 
西 就 会 更 快 落下 。 


刚 开 始 运动 时 以 恒 
定 的 加 速度 9 下 落 ， 
但 是 加 速度 逐渐 变 
小 ， 最 后 物体 以 恒 
定 的 速度 下 落 。 





这 个 恒定 的 速度 由 物 
体 的 重量 决定 。 
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| 与 假设 不 存在 重力 时 一 样 ， 考 虑 初速 度 不 为 0 而 为 w 的 情况 





COREY 
v 


) 代入 v(0)= vo 进行 求解 
CORNE ETIS 
v 

求 未 知 的 常数 


ut)= Lafu- 2) 


v 


= Ze) + we | anaes att 





虽然 从 刚才 的 图 中 可 以 
看 到 ， 若 初始 速度 vo 为 不 管 是 在 哪 种 条 件 下 ， 都 
正 ， 就 会 向 下 移动 ， 初 趋 近 于 终端 速度 mg6rmrl 
始 速度 w 为 负 ， 就 会 向 


气 中 以 初速 度 w 下 落 的 物体 的 速度 随时 间 变化 
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那么 ， 到 现在 为 止 ， e « Š 确实 …… 从 来 也 没有 
落下 来 过 呀 ! 
























































AT EERE ZB m ‘ 
不 下 来 的 原因 了 。 ERATZ REN 




















































































































































































































































































































































































































































































































WTR AREER 用 水 的 密度 1 x 10skg/m? 能 够 计 
能 达到 10hm=(0.01mm)。 算出 终端 速度 。 





110° kg/m? x An(10pm)?x 9.80m/s? 
3 _____ =. 2x10 m/s = 1.2cm/s 
6r x1.8x10°Pa's x 10um 


HL, BREM KN 
以 每 秒 l1cm 左 右 的 速度 
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如 果 上 升 气流 的 速度 大 
于 这 个 速度 ， 


虽然 水 滴 在 空气 中 
下 落 ， 却 整体 向 上 


云 恬 结 核 在 慢 慢 地 下 落 ， 但 是 周围 的 
空气 以 更 快 的 速度 上 升 ， 


所 以 云彩 是 不 会 落下 
来 的 。 





SOA EHRR ME RERA ES 
BR BARD 








就 是 那么 回 事 ! 











到 此 为 止 ， 我 们 以 云 攻 结 核 为 例 ， 
研究 了 只 考虑 粘性 阻力 的 情况 下 小 
物体 的 落体 运动 。 

















两 消 、 球 和 人 等 常 | 
见 的 物体 落下 的 情 


况 是 什么 样 呢 ? 


只 要 考虑 惯性 阻力 
就 足够 了 。 








只 考虑 惯性 阻力 的 情况 下 的 运动 方程 
由 于 重力 与 阻力 达到 平衡 ， 加 速度 dv/dt 
最 终 变 成 0， 只 考虑 惯性 阻力 时 的 终端 速 


) Bv.» 


只 考虑 惯性 阻力 情况 下 的 终端 速度 





可 分 高 变量 分 解 右边 分 数 


ra ea 


(li 
2g\v.-v v +v. 


分 离 变量 进行 积分 


t= Fe (lnle,-osinle,+0))+C 


-mc 
2g |v,-v 
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想 得 知 静止 物体 开始 落下 时 的 样子 ， 设 t = 0 时 刻 时 v(0) = 0， 那 么 积分 常数 
就 变 成 如 下 所 示 。 








由 于 物体 的 速度 v 不 可 能 大 于 终端 速度 v。 ， 所 有 去 掉 绝对 值 符号 也 没有 关 
系 ， 最 后 所 求 得 的 速度 如 下 所 示 。 





RA BREED RAI 


t 
0 


初速 度 为 0 时 ， 空 气 中 下 落 物 体 的 速度 随时 间 的 变化 ( 惯性 阻力 的 情况 ) 





虽然 曲线 与 恭 性 阻力 非常 相似 ， 但 是 这 个 
曲线 更 陡 一 些 。 








tanh x = sinh x/coshr. 
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典型 雨滴 的 终 六 速度 是 
这 样 的 。 
@ 设 半径 =1000 um (=1mm) 
@ 大 气 的 密度 p=1.2Kg/m 
重水 的 密度 p'=1 x 103kg/rm 
y 4 
4x—nx1x10°kg/m’x (1mm)*x 9.80m/s? 
C a Tx1.2kgm°x (Imm)? i 
l | 
CAI 
A 
yr An 
ie 
6.6m/s! 
这 么 说 来 也 没有 





什么 了 不 得 。 













































































































































所 以 就 算 被 两 去 D D 
中 ， 也 不 会 有 生 [lg 
与 云 桥 结核 1.2cm/s 的 如 TERA: 
度 相 比 ， 感 觉 非 常 快 ， 
但 即使 是 以 这 样 的 速度 下 被 十 注 击 不 wR Lo) 
落 400m 也 需要 1 分 钟 以 上 会 多 伤 ! J 
的 时 间 。 
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我 们 着 手 解决 “ 增 减 


FR” HERRED 
模型 。 





| 5. 常数 变易 法 | 


让 我 们 再 回顾 一 下 如 何 推导 出 考虑 黏 性 阻力 时 ， 物 体 落体 运动 微分 方程 的 解 。 


Sag 一 要 求解 的 油分 方程 an 


由 于 不 是 可 分 离 变 量 微分 方程 ， 因 此 无 法 求解 ， 那 么 暂且 忽略 右边 第 一 项 中 的 g 
(暂时 设 为 没有 ) ， 那 么 就 可 以 变 成 可 分 离 变 量 微分 方程 


Bow ATARE (42) 


就 可 以 求解 了 。 得 到 的 解 为 
ulb -ce lf (43) 
由 于 含有 任意 常数 c ( 求 n 阶 微分 方程 的 解 时 ， 不 适用 初始 条 件 等 条 件 并 含有 n 个 任 
意 常数 的 解 叫 做 通 解 ) ， 设 初始 条 件 为 :=0 时 v(0)=v。， 可 得 到 特 解 
v(t)= ne" 1 (44) 
(在 通 解 的 任意 常数 中 代 人 特定 数值 得 到 的 特定 的 解 叫做 特 解 ) 。 若 想 求 可 分 离 变量 
微分 方程 (4.2 ) 的 解 ， 那 么 到 此 为 止 就 可 以 结束 了 ， 但 是 现在 想 求解 的 是 微分 方程 
(4.1 ) ， 需 要 对 可 分 离 变量 微分 方程 的 解 进行 补充 。 假 设 通 解 ( 4.3 ) 的 任意 常数 c 为 
时 间 / 的 函数 c()， 则 函数 化 为 


olt) =t 二 假设 任意 常数 c 为 时 间 / 的 函数 c() (45) 
把 假设 的 解 代入 原来 的 微分 方程 (4.1) 中 就 能 确定 函数 cl), 得 到 解 
的 = 了 + ce- 。 生 想 要 求解 的 微分 方程 的 角 (46) 
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现实 世界 | 数学 世界 


进行 两 次 循环 之 后 得 到 期 待 的 结果 
4 考虑 车 性 阻力 情况 下 物体 在 空气 中 下 落 的 模型 


虽然 已 经 得 到 方程 的 解 ， 但 也 会 对 这 个 解法 的 普遍 性 有 所 顾虑 。 那 么 ， 我 们 试 着 
把 这 个 方法 一 般 化 '。 
用 一 般 函 数 p(z) 和 qlz) 表 达 要 求解 的 方程 (4.1) ， 则 可 以 写成 


D syag) 一 非 同 次 方程 式 (47) 


在 考虑 蔡 性 阻力 时 物体 在 空气 中 的 落体 运动 方程 (4.1 ) 中 ，p(z) 和 q(z) 都 是 常数 ， 
但 是 一 般 来 说 表示 为 z 的 函数 较 好 。 一 方面 ， 没 有 重力 的 情况 下 的 运动 方程 (4.2 ) 变 成 


型 +ptay=0 。 ， 同 次 方程 式 (48) 


虽然 两 个 都 是 线性 微分 方程 ， 但 是 含 不 含 g(z) 项 会 有 所 不 同 。 若 ( 4.8 ) 中 没有 


dz) 项， 微分 方程 中 关于 y 和 -和 的 各 项 具有 相同 次 数 ， 这 种 形式 的 微分 方程 叫做 齐 次 方 
程 。 式 (4.7 ) 中 ,由 于 含有 e(z) 项 ， 不 能 变 成 相同 次 数 ， 因 此 叫做 非 齐 次 微分 方程 ， 
qlz) 叫 做 非 齐 次 项 。 


1 这 也 是 数学 的 优势 之 一 。 


设 积分 常数 为 C， 对 齐 次 方程 (4.8 ) 进行 变量 分 离 ， 可 得 到 
aera e ehed 人 熏 齐 次 方程 的 通 解 (49) 
这 是 含有 任意 常数 的 通 解 。 用 变量 x 的 函数 c(x) 蔡 换 齐 次 方程 的 通 解 中 包含 的 积分 


常数 C， 则 


y= Oho ”+ 用 函数 替换 齐 次 方程 通 解 中 的 常数 
另外 ， HATAR, Meike, M 

y=c(ae eX 全 假设 的 非 齐 次 方程 的 解 (4.10) 
把 上 式 带 入 非 齐 次 方程 (4.7 ) ,求解 假设 的 函数 c(z) 

de(z)e > 


at p(z)c(z)e m= g(x) 





Sez) 
LE tents 2) + nla)eta)e HH = qo) 


SAE) even o(a (plae + pee he qla) 


sa ehad = q(z) 





wo, qirje 
nato 
e(z) aj qtajelneuzdz+e' 人生 假设 的 函数 
把 所 求 得 的 假设 的 函数 cx) 带 入 假设 的 非 齐 次 方程 的 解 (4.10 ) 中 ， 则 
v= (| qerowaz+ ojeirowe tarar a 


这 样 就 求 出 了 此 非 齐 次 方程 的 通 解 。 
那么 ， 展 开 非 齐 次 微分 方程 的 通 解 (4.11 ) ， 则 


y= eli pet | qlaje zdr 十 ciejplzdz 
由 于 c'=ec， 得 知 上 式 为 齐 次 方程 的 通 解 (4.9 ) 


celra +- 齐 次 方程 的 通 解 
与 非 齐 次 方程 的 特 解 
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Epes | daei 。 一 非 齐 次 方程 的 特 解 


的 和 。 也 就 是 说 ， 非 齐 次 方程 的 通 解 是 齐 次 方程 的 通 解 与 非 齐 次 方程 特 解 的 和 。 


4 非 齐 次 方程 的 解 


确认 一 下 。 齐 次 微分 方程 (4.8 ) 的 解 为 y=u(z), 则 


du(z) 


Tas on0 (4.12) 
成 立 。 另 外 ， 若 得 知 y=v(z) 为 非 齐 次 方程 (4.7 ) 的 某 个 解 ， 则 
du(z) (4.13) 





+p(z)v(z) = a(x) 
dz 


成 立 。 分 别 将 齐 次 方程 关系 式 (4.12) 和 非 齐 次 方程 关系 式 (4.13) 的 两 边 相 加 ， 得 到 


dl 
OED ayua) utz) =a) 


因此 得 知 y=u(z)+v(z) 为 非 齐 次 方程 (4.7 ) 的 解 。 通 过 这 个 解法 就 能 够 求 出 非 齐 次 
微分 方程 的 解 。 

这 种 解 非 齐 次 方程 的 方法 叫做 常数 变易 法 。 为 求 非 齐 次 微分 方程 的 解 ， 将 齐 次 
方程 通 解 中 的 常数 进行 变易 ( 把 非 齐 次 方程 的 解 写 成 齐 次 方程 的 通 解 x 未 知 函 数 的 形 
式 ) ， 由 此 产生 了 这 样 的 叫 法 。 





© 求解 非 齐 次 方程 的 方法 


运用 常数 变易 法 时 ， 首 先 在 要 求解 的 微分 方程 中 ， 求 出 求解 较 容易 部 分 的 微分 方 
程 ( 齐 次 方程 ) 的 解 ， 之 后 对 其 进行 补充 ， 以 满足 要 求解 的 微分 方程 ( 非 齐 次 方程 ) 
的 求解 方法 。 先 不 要 被 微分 方程 的 形式 所 限制 ， 从 较 容易 求解 ( 可 求解 ) 的 部 分 开始 
求解 ， 后 续 对 细小 的 部 分 进行 修正 即 可 。 这 是 一 个 非常 好 用 的 方法 。 比 如 空气 中 物体 
落体 运动 的 例子 ， 随 着 时 间 的 变化 物体 以 指数 函数 的 方式 运动 ， 因 此 要 先 求 出 这 个 部 
分 的 解 ， 然 后 再 进行 补充 以 满足 要 求解 的 微分 方程 ， 让 一 切 都 合乎 轴 辑 即 可 。 
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1 振动 现象 












































“a 
参拜 几 次 愿望 就 
会 实现 ! 
这 个 神社 真是 
太 棒 了 ! 
N Z 
= am, \Y In 5 
= \e Gu) trenk 
一 faite À 
fi 
这 个 家 伙 


RECTED! Rey 


on 








来 这 里 参天 的 
话 ， 总 觉得 会 
变 得 更 加 精力 
充沛 啊 。 



































































































































啊 ， 明 白 了 。 

















lol es 
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` Mea 
vs moai 
evaj i 





虽然 外 观 不 太 
好 看 





得 WY, 





raans, se W 
Fl 得 会 很 有 帮助 的 。 
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一 会 儿 增 加 一 会 儿 减 少 
EBB AF 


振动 现象 需要 用 比 以 往 更 加 复 
oe “二 阶 线性 微分 方程 ”来 


这 最 后 的 内 容 ， 大 家 
MEM SHA AA! 


BNA KZN 
HAE 


骑 自 行车 的 时 候 屁股 
不 会 疼 也 是 因为 弹 黄 


弹 黄 装置 通过 反弹 
来 续 和 冲击 。 








| omman | vanes, 


那么 ， 这 个 装置 用 


FER HK 
长 度 。 











KAHL? 














不 附加 任何 外 力 状态 下 的 
用 手 找 一 下 小 球 再 松手 ， 
弹 黄 就 开始 振动 。 弹 黄 长 度 叫 做 自然 长 度 。 



































































































光滑 的 水 平台 










































一 一 一 
a Cee (ae 
BRAT, aL 
不 回去 了 。 
"ey 但 是 在 一 定 范围 内 ， 
| aD 拉力 和 收缩 力 保持 平 
衡 时 ， 就 能 恢 能 复 到 
自然 长 度 。 
看 看 小 球 反动 的 样子 
Bhatth 2 E Q exe 
BT! 
mn 弹 英 的 伸 长 量 与 拉 
i 力 之 间 有 什么 样 的 
* ç 2 关系 呢 ? 
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收缩 、 拉 仲 说 起 来 怪 
麻烦 的 ， 要 不 把 收缩 
叫做 “反方 向 拉 伸 ” 
吧 


那么 ， 弹 黄 保持 自然 长 度 时 ， 小 球 
的 位 置 叫做 原点 ， 设 小 球 的 伸 长 量 
为 Y， 则 可 以 得 到 如 下 的 关系 式 。 















































照 你 这 么 说 


弹 鞭 的 收缩 力 与 收 
缩 量 也 是 同样 的 关 


这 个 关系 式 叫 做 胡 克 
定律 ， 所 受 的 力 叫做 
弹力 。 





小 球 的 受 力 方向 是 拉 伸 的 相反 方 
页 








因 弹 黄种 类 而 不 同 ， 称 为 
弹性 系数 。 












































































































































































































































X 值 大 的 弹 黄 和 小 的 弹 黄 进 行 
同样 拉 伸 时 ， 


RK OHH ER ARK 
的 拉力 。 


也 就 是 说 ， 弹 性 系数 表 
THER, 


作用 于 弹 黄 的 力 不 
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阻力 以 空气 阻力 和 摩擦 力 
的 形式 作用 于 运动 的 相反 
方向 。 


这 里 ， 
BMRA BSR 
速度 成 正比 的 阻力 。 


车 用 扩 表 示 阻 力 ， 就 可 以 用 小 
球 的 位 置 随时 间 的 变化 表示 小 
球 的 滑动 速度 。 
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[Raho Toe - 


如 果 不 持续 施加 外 力 
不 要 拉 一 次 弹 黄 就 放手 ， 的 话 ， 弹 黄 一 会 儿 就 
需要 持续 施加 外 力 。 




















过 一 会 儿 弹 芙 就 被 拉 坏 了 ， 
所 以 一 般 来 说 要 施加 周期 性 
的 外 力 。 

















暂且 把 它 看 做 
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正如 目前 为 止 所 看 到 
的 ， 这 个 振动 体系 ， 









不 多 外 力 的 情况 下 ， 
物体 保持 原来 运动 
状态 性 质 。 


30 
动 





TARRE, RR, Mh 


三 个 方面 进行 考虑 。 
































电路 中 电流 的 振动 
也 能 够 用 同样 的 力学 模型 表示 。 











该 模型 不 仅 限于 力学 上 的 振动， 
还 能 表示 更 广泛 的 振动 现象 。 





这 种 情况 下 ， 电 感 (L) 、 
eg (C), 电阻 (R) 等 
电路 的 要 素 ， 


与 力学 中 的 惯性 、 弹 性 和 阻 
力 等 因素 相对 应 。 


那么 ， 写 一 下 这 个 振动 
体系 的 运动 方程 吧 。 


设 小 球 的 质量 为 mm， 那么 方程 | 
可 以 这 样 表示 。 A 






































ae 
mh | BA Hn 
g raat ees 
i ep eee BO) 
| dt? dt | 
3 的 风量 ELIA 
gamhna 














7 





这 样 就 完成 了 振动 
体系 的 建 模 。 
zá 
g y 


那么 ， 解 一 下 这 个 
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[e 振动 模型 2 eia ] 






























































这 个 方程 到 底 从 哪儿 
分 解 好 呢 








比如 刚才 这 种 情况 
RELA? 


m—+c—+ker= F(t) 
d? dt 7 

































































试 着 考虑 一 下 没有 MRK ERE R 
阻力 和 外 力 的 情况 。 性 这 两 个 要 素 咽 ? 








这 么 一 来 ， 只 有 弹 呐 的 
弹力 作用 于 小 球 。 











为 了 使 上 式 看 起 来 
简单 明了 ， 





左边 只 留 一 个 x 的 二 阶 导 数 ， 
将 kr 移 项 。 


就 变 成 这 样 ! 





这 就 是 二 阶 齐 次 
MDF ALE! 











































































































© 整个 流程 就 是 这 样 的 。 刚 开始 考虑 的 模型 较 难 ， 不 知道 如 何 用 数学 知识 处 理 才 
好 ， 通 过 对 现象 进行 简化 、 对 方程 进行 降 阶 ， 就 可 以 重新 考虑 这 个 模型 啦 。 





1 
Li 
¢ 





MREY SUSE 


1 
1 
1 
1 
1 
有 
¢ 








现实 世界 ; 数学 世界 


4 在 没有 阻力 和 外 力 的 情况 下 ， 挫 在 弹簧 上 的 小 球 的 运动 模型 


是 。 
那么 ， 试 着 求 一 下 不 存在 阻力 与 外 力 情况 下 的 通 解 。 请 仔细 观察 这 个 微分 方程 。 
又 看 啊 ? 


要 目不转睛 地 看 ! x 对 + 进行 两 次 微分 ， 虽 说 带 bm 这 个 系数 ， 但 是 再 回 到 -x*， 注 意 
这 个 式 子 !'。 
微分 两 次 就 能 变 回 原 函数 …… 好 像 在 哪里 看 见 过 呀 ! 


O @ese 


1 不 认真 看 ， 是 不 会 注意 到 的 。 
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O STERKAR 
A i =cost 
dt 


再 一 次 对 /进行 微分 ， 再 套用 余弦 函数 的 微分 公式 


i cost=-sint 


得 到 
2 
T i = A cost =-sint 
dé? dt 


这 样 两 次 求 导 之 后 变 成 原 函数 。 巧 妙 运用 三 角 函 数 的 这 种 性 质 ， 就 能 够 得 到 方程 
的 解 。 不 过 ， 有 人 会 考虑 在 二 阶 齐 次 微分 方程- 秘 -=- 去 z 中 ， 有 klm 这 个 系数 ， 
该 怎么 处 理 呢 。 利 用 微分 的 基本 性 质 之 一 ， 设 为 常数 ， 则 

过 ints coats -wzsinwt 

dt? dt 


实际 上 , 设 
Eia 61) 
m 

便 得 知 
a(t) = sinwt (5.2) 


为 微分 方程 4- Le 的 解 。 太 好 了 ! KET! 
不 过 ， 且慢。 就 只 有 这 一 个 解 吗 ? 解 (5.2) 的 整数 信 
z(t) = Asinwt (5.3) 
也 是 微分 方程 =- 直 z 的 解 。 
把 这 个 整数 倍 的 解 代入 到 微分 方程 中 就 能 确认 是 不 是 方程 的 解 。 





EAsinwt = -Phsinot 
ae 


BAN, HERRERIA 
£ cost =-sint 
dt 
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对 4 进行 微分 ， 利 用 正弦 函数 的 微分 公式 

sint = cost 
可 得 到 

Riera Sais =-cost 

dt? dt 
设 B 为 常数 ， 

a(t) = Beoswt (5.4) 
thane Se -Erhi 

实际 上 ， 通 解 就 是 解 ( 5.3 ) 和 (5.4) 的 线性 组 合 ?。 

z(t) = Asinwt + Bcoswt (5.5) 
就 是 没有 阻力 与 外 力 情况 下 的 通 解 。 
将 通 解 (5.5 ) 对 :进行 两 次 微分 


2 
E (sinut + Bcoswt) = -w*(Asinwt + Bcoswt) 


就 可 以 确认 是 微分 方程 的 解 。 


2 常数 倍 求 和 得 到 的 。 











RRA TOM, 
那么 通过 研究 具体 的 状况 
解释 一 下 这 个 解 吧 。 





没有 限 力 与 外 力 情况 下 的 擅 解 





a(t) = Asinwt + Bcoswt 





HORA BAF. 





拉 住 小 球 把 弹 黄 伸 长 至 To 的 
位 置 。 























初始 条 件 变 成 (0)=Zo， 
[dz/dt]o=0 





z(0)= Asin(w-0) + Bcos(w-0) = B 





[canine + Bcoswt) | =wAcos(w0) - wBcos(w:0) =wA 


= 














z(t) = x coswt 
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这 个 解 表示 小 球 从 -zo 
的 位 置 到 zo 的 范围 内 周 
期 性 振动 的 状态 。 





是 单纯 的 振动 呀 。 


周期 是 振动 一 次 所 需要 
的 时 间 ， 我 们 设 振动 的 
AMAT. 








意思 是 说 ， 
解 中 coswt 的 wi 要 成 为 2x 时 ， 





也 就 是 说 ， 若 没有 阻力 个 _ 


和 外 力 ， 会 一 直 保持 国 直 
RBE g 













































































像 这 样 振幅 和 周期 
都 国定 的 振动 ， 


另外 ，w= kim 
叫做 这 个 振动 体系 
HOA ARR. 








BEDS 
HESS, 
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© 现在 的 情况 如 下 所 示 。 


MBS SU as 











现实 世界 数学 世界 
通过 降 阶 能 够 得 到 简化 之 后 的 模型 的 解 。 现 实 中 ， 阻 力 必然 存在 ， 所 以 说 的 是 阻力 非常 小 的 情况 。 


没有 阻力 与 外 力 的 情况 下 ， 挫 在 弹 赞 上 的 小 球 的 运动 


简 谐振 动能 够 解释 这 里 所 提 及 的 忽略 阻力 情况 下 弹簧 振动 和 振幅 小 的 振子 运动 。 
另外 ， 了 解 这 些 概念 是 非常 重要 的 ， 因 为 它们 是 研究 振动 现象 的 基本 知识 。 
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L 4 振动 模型 9 有 阻力 的 情况 ] 
有 阻力 时 的 解 














那么 ， 若 考虑 阻力 ， 
情况 会 是 什么 样子 
呢 ? 








两 边 除 以 m 进 行 整理 ， 
则 变 成 这 样 。 
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这 就 是 到 现在 为 止 的 流程 。 接 下 来 我 们 把 降 阶 简化 之 后 的 二 次 循环 模型 上 再 加 上 
别 的 因素 ， 使 其 变 得 稍微 复杂 一 点 。 







SBMRSYSRUBAIY 





现实 世界 数学 世界 


4 有 阻力 时 ， 挫 在 弹 笑 上 的 小 球 运动 的 模型 化 


© 那么 ， 求 不 存在 外 力 时 的 通 解 。 这 时 也 要 注意 观察 这 个 微分 方程 


Le + 24, wr=0 


© 好 的 ! 


© 那么 ,我 们 能 够 注意 到 ，z 对 4 进行 两 次 微分 得 到 的 项 、 进 行 一 次 微分 得 到 的 项 与 = 
本 身 相 加 为 0。 也 就 是 说 ， 这 里 的 函数 z(t) 必 须 是 微分 后 形式 不 发 生 改变 的 函数 。 





3 由 于 在 简 谐 运动 的 时 候 锻炼 了 眼睛 ,应 该 比 以 前 更 早 看 到 了 RE 
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微分 之 后 形式 不 发 生 改 变 的 函数 ， 似 乎 在 什么 地 方 见 过 。 
这 也 是 在 53 页 出 现 过 的 式 子 ， 观 察 指数 函数 的 微分 公式 


Aa 
dt 


之 后 就 能 明白 ， 指 数 函数 就 是 进行 微分 之 后 形式 不 发 生 改 变 的 函数 。 巧 妙 地 利用 
指数 函数 的 这 个 性 质 ， 能 否 求 得 这 个 微分 方程 45 42,22 +u2z =0 的 解 呢 ? R 
们 试 试看 吧 。 

指数 函数 的 微分 公式 中 ，e 对 ;进行 微分 之 后 跟 以 前 一 样 ， 也 不 会 出 现 系数 什么 
的 。 微 分 多 少 次 都 是 相同 的 形式 。 再 看 看 这 个 微分 方程 -4 +2722 4 ute 0, 
式 子 中 各 项 都 含有 系数 ， 就 这 样 求解 肯定 不 可 以 ， 因 此 需要 动 动脑 筋 ， 设 为 党 
数 ， 结 合 微分 的 基本 性 质 (53 页 )* 得 到 


Sew =Ae™ 
dt 





通过 反复 微分 可 以 得 到 
Zo =r 


出 现 了 系数 ， 可 以 用 于 微分 方程 。 


这 里 ， 我 们 尝试 这 样 一 种 方法 ， 把 想 要 求解 的 微分 方程 43 + 2 和 + wzz=0 的 解 
假设 为 

z(t) =e" “< 假设 的 解 (5.6) 
WRB A Ta Ee $27 Shute =0 方程 的 话 ， 确 定常 数 入 。 


dt 
把 假设 的 解 (5.6 ) 代入 微分 方程 =. +i ute =0 


CRETE AE EA 
— e + 2y — e tuen 0 (5.7) 
dt? dt 
再 进行 微分 计算 

Me™ + 2ye” + wzex=0 
合并 同类 项 ex 

(2+ 27A +u*)e*= 0 








4 这 叫 复合 函数 的 微分 。 


由 于 ex 不 可 能 为 0， 为 满足 式 (5.7) ， 则 必须 有 


+ 27A+u7=0 (5.8) 


这 是 关于 的 二 次 方程 啊 。 
那么 ， 这 个 微分 方程 的 解 为 


=F t-a, eee ar (5.9) 


这 就 是 说 ， 与 代数 方程 (5.8 ) 的 两 个 解 (5.9 ) 相对 应 ， 可 以 很 巧妙 地 求 得 172 页 


OO 


的 微分 方程 .4 + 24 SE + ut = OAAR 

















Zi=e rme 全 微分 方程 的 解 (5.10) 
© 但 是 ,这 些 解 都 不 全 任意 常 数 ， 很 站 二 都 不 是 通 解 。 但 是 这 时 我 们 知道 应 该 怎样 
才能 解决 问题 。 
ORELLA 
© th. 那么 为 求 得 通 解 ， 将 解 (5.10 ) RAER O 
amc (Ne 二 假设 系数 的 解 (5.11) 
那么 ， 使 用 常数 变易 法 ， 假 设 系数 的 解 (5.11) 对 /进行 微分 
drte) _ de, O a, de^ 
dt dt OF 
“ty ) ete, (Ae 
doit) > de (t) g de,(t) des ae sO secon ges 
dt? dt? dt dt 
š Ho, aa ca dapat 


把 上 式 代入 微 分 方程 -4 + 27 SE + ute = 0 中 


5 曲线 与 x 轴 不 相交 。 
6 二 次 方程 az?+bz+c=0 的 求解 公式 为 == (-b2 VF dae )/ 2a。 
7 注意 到 了 吧 ? 这 就 是 忧 然 大 悟 之 处 。 了 前 面 加 上 系数 2， 使 二 次 方程 的 解 更 简单 。 
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(oo, Hase) 
dt? 


s(t ew +e (t)A e> ) +u7e,(the'=0 
再 合并 相同 项 ex 
pau. +20,+ SO. er, oh (5.12) 


(4) 不 要 被 外 观 的 复杂 所 困惑 ， 这 里 的 e“ 同 样 也 不 可 能 为 %， 为 满足 式 (5.12) 
则 


#20 ao DE (tm) -0 


另外 ,根据 (5.8) ， 有 +27A+ w=0， 则 


de(t) doi(t) 
dé 


20, + — =0 
+20,+ 9): dt 


PAERD, es 


enne 一 < Be, D 20,4 gean de, so =0 


“en 


这 样 ， inner 另外 ,根据 代数 方程 的 解 (5.9 ) 


Aart Pw, yen Pw 


求 和 入 ,的 和 
At dye (Fo) 
利用 这 个 关系 ， 则 


(er ) =0 成 立 。 
设 C 为 任意 常数 ， 然 后 将 上 式 对 :进行 积分 ， 


8 这 种 方法 有 点 技巧 性 ， 但 较 常用 。 
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Se 


erat 


由 于 dag 
at 


设 积分 常数 为 c， 想 要 求 的 函数 c(f) 变 成 


a(t) -cf Oat +a (5.13) 
那么 ， 这 个 积分 必须 分 为 A,= 和 和 和 ,和 两 种 情况 考虑 ， 随 后 分 别 进行 讨论 。 
首先 ， 考 虑 一下 NM- 的 情况 。 这 种 情况 下 ， 由 于 eorsx=1， 函 数 (5.13) 变 成 

ot9-c| dt+a 
积分 运算 ， 则 

aO-Ct+ra 


把 现在 求 得 的 函数 c,( 1) 代 和 人 假设 未 知 系数 的 解 ( 5.11 ) 中 ， 就 能 够 得 到 微分 方 
BEE 42 awe 0 的 通 解 


a(t) = (Ct + oa)e™ 


根据 代数 方程 的 解 ( 5.9 ) 

Am-yty P-a, NNY- 
入 = 和 ， 也 就 是 说 

yt Pew =y- Pw? 

YY-w =0 
we? 

这 时 由 于 

N-X--7 


把 常数 C 换 写成 6， 结 果 所 求 得 的 通 解 为 
z(t) = (a + Bt)e™ uP FR F AE 
这 里 ，a 和 [8 均 为 任意 的 常数 。 
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© 接 下 来 ,考虑 A, 关 和 的 情况 时 ， 由 于 函数 ( 5.13 ) 为 指数 函数 的 积分 ， 进 行 积分 运 
算 , 则 


e(0-—S 


和 -入 


这 里 ， 代 数 方程 的 解 (5.9 ) 为 
dant P-a, M=- VP - 


DAFA LIRA BIT UE Bla? é 7, ARS PE SP 
wo, we? 
两 种 情况 进行 考虑 。 


COM a 





@ 首先 考虑 -下 wz> ?的 情况 。 这 种 情况 下 ， 代 数 方程 的 解 (5.9) 
Amy y P-a, NNT 
的 根 号 中 的 值 会 变 成 负数 ， 此 时 的 解 为 复数 。 引 入 虚数 单位 ;, 则 可 写成 


M+iQ, A,=-7-i2 


anyag (5.14) 
设 a 和 /为 任意 的 常数 ， 因 此 所 求 的 通 解 为 
X(t) = aei + Bei 一 必 > 了 情况 下 的 通 解 


那么 ， 必 < 他 的 情况 会 是 什么 样 呢 。 这 种 情况 下 ， 代 数 方程 的 解 ( 5.9 ) 为 
Amy y P-a, NT 
变 成 实数 ， 则 写成 


à= -y+T, N=-7-T 
ra{7-w 


设 a 和 /为 任意 的 常数 ， 因 此 ， 所 求 的 通 解 为 
a(t) = aen 4 Beer oF < PRO FB 





ETETE? i 






























































































































wF> 7: a(t) = ae M4 Be, Quy u- 





考虑 阻力 的 情况 下 ， 
微分 方程 的 通 解 可 分 











w=: c(t) = (a + Bt)e™ 
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图 考虑 阻力 情况 下 的 解 1 (阻尼 振动 ) 
| 首先 ， 研 究 一 下 wW2 > 人 的 情况 








O 设 < 和 为 任意 的 常数 ， 那么 ， 这 种 情况 下 的 通 解 为 


a(t) = ae + eto, Qayu- 


将 指数 部 分 分 为 实数 部 分 和 虚数 部 分 ， 则 
a(t) = e ”(ae' + Bei) (5.15) 


© 那么 ， 通 解 (5.15) 的 括号 中 的 aeie+ pet, HARRA 
etz= cosZ 土 isin7 
则 
a(t) =e {a(cos Qt + isin Q) + Blcos Qt - isin 2t)} 
= e{(a + )cos Qt + i(a - Asin Qt} 
=e (acos Qt + bsin Qt) 
{Bika=atf, b=i(a-). 
®© 虽然 此 解 与 周期 为 2w/w 的 简 谐 振动 的 通 解 (5.5 ) 非常 相似 ,但 是 有 两 个 不 同 
点 ,一 是 周期 为 2n/Q 和， 二 是 整体 有 e-* 这 个 指数 函数 相 乘 。 这 个 解 所 要 表现 的 
是 什么 样 的 运动 呢 ? 
与 简 谐 运动 时 一 样 ， 拉 住 小 球 把 弹簧 拉 伸 至 ro，t=0 时 开始 放手 。 假 设 初始 条 件 为 
t=0 时 ，z(0)=z。，[dz/dt],-。=0, 根 据 式 (5.16) 则 


(5.16) 


z(0)=e {acos(20) + bsin(2-0)} =a =z, (5.17) 


9 实际 上 ， 这 是 将 指数 函数 扩大 到 复数 情况 下 的 定义 ， 因 此 反 过 来 说 也 是 成 立 。 


a=zu， 将 式 子 (5.16 ) 对 t 进 行 微分 ， 由 于 a=z。， 则 


£ = ye" (acos Qt + bsin Qt) + e(-aQsin Qt + bQcos Qt) 
= e™{(-ya + bQ)cos Qt - (ya + aQ)sin Qt} 


根据 上 式 得 到 
加 -emof(-Ya+bDjcos(20)-(Oa+a9jksin(20)} 


=-ya+bQ=0 (5.18) 
abelar, 


Q 
最 后 ， 根 据 式 (5.17) 和 (5.18) 得 到 解 


a(t)= enz, (cos ot + Asin 21) . 
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从 图 上 看 出 ， 小 球 以 
27/0 为 周期 振动 的 
同时 ,振幅 也 逐渐 地 


让 我 们 看 看 连接 扰动 的 
各 顶点 的 线 。 








ARKAT. 





a <y RR FAR (BERN) 














变 成 指数 函数 ec， 






VRRP, 
决定 振幅 减 小 的 速 
度 。 





Jit RRR 周期 与 相对 应 的 时 间 常 数 之 
数 T=1/7 叫 时 比 TI7 = 20y/Q, ATR 
间 常数 。 RAF, 











HPA 








TIT =2ry/Q 








阻尼 振动 的 振幅 就 变 
BARR HTH, 











阐 谐 振动 的 周 





有 阻力 情况 下 的 
振动 比 简 谐 振动 
稍微 继 慢 一 些 。 
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那 我 们 研究 一 下 自行 车 弹 黄 装置 中 的 弹 黄 
阻尼 振动 的 实例 。 总 是 在 振动 着 ， 





























那么 ， 在 自行 车 上 
安装 一 个 叫做 减 震 
器 的 调整 阻力 的 装 
置 ， 就 能 达到 快速 






超市 或 者 便利 店 里 的 动 
门 上 也 在 使 用 减 震 器 。 





































































































如 果 没有 减 震 器 Š 
是 吧 哄 吧 哄 的 响 ， 很 
危险 的 呀 。 
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图 有 阻力 情况 下 的 解 2 ( 过 度 衰减 ) 
AL 接 下 来 ， 我 们 研究 一 下 ?< 的 情况 


合 。 设 a 和 为 任意 的 常数 ， 这 种 情况 下 的 通 解 为 
z(t) maen + Ber”, T= y P-a 








(5.19) 

那么 像 往常 一 样 把 指数 部 分 分 为 两 个 部 分 ， 则 
a(t) =e"(ae"+ Be”) (5.20) 
从 这 个 式 子 并 不 能 判断 出 哪 一 项 与 振动 有 关系 %。 那 么 ， 此 解 所 描述 的 是 什么 样 


的 运动 呢 ? 


®© 拉 住 小 球 把 弹簧 拉 伸 至 z。，t=0 时 开始 放手 。 初 始 条 件 为 {=0 时 ，z(0)=z。，[dz/ 
dt) -=0， 根 据 通 解 (5.20) ， 则 
20) =e (ae Be) =a+ B=x, (5.21) 
另外 ， 通 解 (5.20) 对 + 进行 微分 
a = -entaen+ Ber”) +eM(ale”- Bre-n) 
=e{('-yae*- (I+ 7) he} 


根据 上 式 ， 则 化 为 
dz) aeri- yae- (r. -roy-0 
Ge ET- Wae™-(T+ 7) Be} ua 
“(T-Wa-(T+W8=0 
最 后 ,根据 式 (5.21) 和 (5.22), W 


1 不 管 是 增加 的 还 是 减少 的 ， 都 是 单调 函数 。 为 了 表达 振动 的 函数 ， 指 数 部 分 必须 要 引入 


第 5 章 %9 二 阶 线 性 微分 方程 185 


a+B=z0 
T-a (P+ B=0 


saghe 


因此 ， 解 为 
omer [Bheir] 
虽然 这 样 也 可 以 ,但 是 用 双 曲 线 函 数 " 
Ete Cag? 


, sinhz= 
2 


的 话 ， 可 以 得 到 与 欧 拉 公式 相似 的 关系 
er=coshz+sinhz 
利用 上 式 ， 则 通 解 (5.20) 为 


coshz= 








a(t)= fih 7) (cosh Tt+sinh Ft) + 2i- 7) (cosh Tt- sai) 


=e%z, (cosh Ti+ Tsinh r) 


这 样 就 变 得 稍微 简单 些 了 。 


11 读 作 正弦 双 曲 线 和 余弦 双 曲 线 。 
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小 球 的 位 移 不 可 能 
Too 


WP<J WR FAR (HARB) 








a ar, 起 建立 扰动 体系 的 楼 
< 型 ， 才 发 现 原来 也 有 


不 进行 振动 而 发 生 
表 碱 ， 这 样 的 运动 
AERA. 


制作 模型 的 时 候 ， 
即使 是 没有 设想 过 的 
RE, 


通过 微分 方程 的 形式 引入 数 
学 世界 时 ， 也 能 够 包含 所 有 
的 情况 。 


这 也 是 利用 数学 的 优点 
之 一 ! 
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过 度 表 碱 的 情况 通 解 是 两 个 指数 函数 


F, RABE 的 和 RET. 
度 取决 于 7 的 大 小 。 
PR 








a(t)= ae petn, 了 = 他- 二 

























Faye 


aoe ; 
Ee 
1 变 大 ， 则 立即 发 生 表 减 
NK 


那么 ， 根 据 这 个 效果 ， 
RABE KA Ra, 





根本 不 发 生 表 减 。 




















自行 村 的 缕 冲 装置 
调整 过 度 表 减 时 ， 


如 果 人 突然 离开 车 
È, FÈRRE 
振动 ， DRAB 
升 起 。 


Fay RRM 
时 ， 需 要 的 时 间 是 择 
动 门 完全 关上 为 止 的 
时 间 。 


这 种 运动 好 像 也 
不 能 运用 到 什么 
Arn? 





沉 下 去 之 后 ， 
FEBRRE. 


可 以 用 于 高 级 产品 。 


也 许可 以 面向 高 级 音响 设 
备 ， 如 CD 和 DVD 的 开关 





总 感觉 反应 迟钝， 


那么 空调 的 效率 会 
变 得 非常 差 的 
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图 有 阻力 情况 下 的 解 3 ( 临界 衰减 ) 
| 这 次 ， 考 虑 w2= 7 的 情况 








合议 c 和 6 为 任意 的 常数 ， 这 种 情况 下 的 通 解 为 
a(t) = (a + Bt)e* (5.23) 
在 这 种 情况 下 ， 也 不 能 判断 出 与 振动 有 关系 的 项 。 
©) 与 之 前 的 情 吕 一样， 拉 住 小 球 把 弹 和 拉 伸 至 z，，t=0 时 开始 放手 。 初 始 条 件 为 t=0 
时 ，z(0)=z。，[dz/dt],-。=0， 根 据 通 解 (5.23) ， 则 
z(0)=(a+ B-0)e=a=2, (5.24) 


由 于 上 式 ,a=z,， 通 解 (5.23 ) 对 /进行 微分 
Fw Bes (a+ BE 
= {8-a + Aoje” 
根据 上 式 ， 化 为 
图 =(8-(a+80)}e7=0 
“ibe, (525) 
那么 根据 (5.24)、(5.25) 解 为 


z(t) = (zo+ 120t)e™ 


sai (5.26) 
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这 种 情况 下 ， 小 球 的 位 
置 也 不 可 能 超越 原点 ， 











TAL FF a RM, 









































如 果 只 看 临界 表 减 
曲线 的 形状 ， 与 过 
度 表 减 并 没有 什么 
Bal, 














这 里 的 临界 表示 
边界 。 








但 是 临界 表 减 有 更 重 
要 的 意义 。 
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若 阻 力 比 临界 表 减 
小 ,就 变 成 阻尼 振 
动 ， 阻 尼 振 动 的 情 
ATF, 





明白 了 吗 ? 



















































































































































































































































































Shit RAGA MRR 
ARRAN. 














ap 
9 


AY 








有 一 个 叫 自动 闭 门 器 的 
生活 用 品 。 














住宅 的 门 上 也 
经 常 使 用 这 个 
未 西 。 






MERA RE 
器 使 其 接近 临界 表 减 。 






ERRUZ 
RIRTAB. 














如 果 将 其 调整 为 阻尼 相反 ， 
拔 动 ， 





FALARA, 
NARRE 


AB 2 HFA ak Ek 
度 最 快 ， 也 许 门 全 
坏 掉 。 


如 果 没有 门框 ， 
没准 门 会 飞 出 去 ! 
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FABRA, ABH 
振动 周期 一 丰 的 时 间 ， 


但 不 会 超越 原点 。 门 会 很 平稳 地 关上 。 


小 球 就 能 到 达 原 点 的 
位 置 。 














: 如 有 果 把 自行 车 的 缕 冲 效 置 
将 计量 表 调 整 为 临界 |] 、 也 调整 成 临界 来 减 ， 会 发 
kà, gataz AN : 挥 较 好 的 效果 。 










































































Shy Bs ae 








Maca, 
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© 今天 就 到 这 里 。 对 没有 外 力 情况 下 的 三 种 解 进行 了 解释 。 
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L sa—sane | 








®© 到 现在 为 止 ， 以 推 在 弹簧 上 小 球 的 运动 为 例 ， 学 习 了 常 系数 二 阶 线性 微分 方程 的 
解法 ， 在 这 里 进行 一 下 总 结 。 为 求解 这 个 常 系数 二 阶 奇 次 线性 微分 方程 


a Wynd REC BNA (521) 
dz? dz 

设 和 为 常数 ， 假 设 解 为 
gz)=ex ”人 一 假设 的 解 (5.28 ) 


把 假设 的 解 ( 5.28 ) 代入 微分 方程 (5.27) 
g d 

ria be fee 
进行 微分 计算 

ade + be + ce¥ = 0 
合并 同类 项 

(aX*+ b+ cle = 0 (529) 
与 前 面 的 情况 相同 ， 指 数 函数 e~“ 在 z 所 有 的 范围 内 都 不 可 能 为 %， 若 要 (5.29) 成 
X, 必须 满足 

aX+bA+c=0 
由 于 这 是 关于 入 的 二 次 方程 ， 那 么 ， 解 为 


he be 4ac 4i Ayt -4ac 
1 > 


也 就 是 说 ， 能 名 得 到 与 代数 方程 (530 的 后 角 相对 上 的 、 微 分 方程 (5.27) 的 两 
个 解 


(5.30) 





yee, Y=e 二 常 系数 二 阶 奇 次 微分 方程 的 解 
像 这 样 通过 代数 方程 的 解 求 得 微分 方程 的 解 ， 那 么 这 个 代数 方程 (5.30 ) 就 叫做 
特征 方程 。 如 果 已 知 特征 方程 的 解 ， 也 就 能 够 得 到 相应 的 微分 方程 的 解 。 即 如 果 
能 够 求解 二 次 方程 ， 也 必然 能 够 求解 相应 的 常 系数 二 阶 奇 次 微分 方程 。 
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©) 根据 根 号 中 判别 式 b>-4ac 的 符号 /一 ( 根 号 )， 把 特征 方程 (5.30 ) 的 解 分 别 为 几 种 
情况 。 分 别 是 ， 太 >4ac 的 情况 下 有 两 个 实数 解 ，b-4ac 的 情况 下 为 重 解 ， 如 >4ac 
的 情况 下 有 两 个 互 为 共 令 的 复数 解 ， 下 面 我 们 分 别 讨论 一 下 不 同 的 情况 。 
刀 >4ac 的 情况 下 ， 由 于 特征 方程 (5.30 ) 有 两 个 实数 角 
4-34, 起 一 - 了 
a 
设 a 和 为 任意 的 常数 ， 那 么 微分 方程 (5.27) 的 通 解 为 


ters 
加 


y(z)=ae a a + Be ™* 


b=4ac 的 情况 下 ， 特 征 方程 (5.30) 的 解 为 重 解 
b 
各 


设 a 和 A 为 任意 的 常数 ， 那 么 微分 方程 (5.27 ) 的 通 解 为 





Wo)-(a+Bre ™ 
如 >4ac 的 情况 下 ,特征 方程 (5. ny APA ASEH NR 
a= aes ate Le PT) 
KIRE Fa (5.27) 的 通 解 为 


bei 


b 
tiu -， 
dj-a" * +fe™ ™, 
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| 6. 再 回 到 振动 模型 1 有 外 力 的 情况 | 








有 外 力 时 的 解 


那么 ， 





考虑 一 下 外 力 。 



































终于 可 以 计算 最 初 提出 
的 微分 方程 了 。 


如 果 以 恒定 的 外 力 
F(D) 作 用 于 小 球 ， 


弹 黄 向 外 力 的 方向 
HAWK, RH 





ARHED 
是 啊 ， 这 就 是 包含 弹 


力 、 阻 力 与 外 力 的 弹 
黄 的 运动 方程 。 





上 述 这 种 情况 并 不 理想 ， 


若 以 周期 性 的 力作 为 外 力 ， 


F(t) =F cosvt 
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© 这 是 包含 弹力 、 阻 力 和 外 力 的 弹簧 的 运动 方程 ， 两 边 除 以 m， 进 行 整理 


d? mdt m m 


© 这 也 是 一 阶 和 线 性 微分 方程 呀 。 
© sam 
Euih 
m 
另外 ,与 172 页 相同 ， 设 
Lay 
m 
5 
fas 
m 
则 微分 方程 (5.31 ) 可 改写 为 
oo oy dt 全 考虑 外 力 时 的 运动 方程 (532) 
rT a 考虑 HANZ: i 
这 就 是 非 齐 次 微分 方程 。 那 么 ， 求 一 下 这 个 微分 方程 的 解 。 


® 由 于 是 非 齐 次 微分 方程 ， 就 像 在 第 4 章 见 过 的 那样 ， 用 常数 变易 法 进行 求解 。 根 据 
172 页 的 内 容 ， 与 非 齐 次 方程 ( 5.32 ) 相对 应 的 奇 次 方程 为 


效仿 180 页 中 心 > 的 情况 下 奇 次 方程 的 解 ， 设 非 齐 次 方程 (5.32 ) 的 解 为 


Z(t)= Acosut+ Bsinvt ”一 假设 的 解 (5.33) 
4 、B 为 满足 微分 方程 (5.32 ) 的 常数 。 
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© 将 假设 系数 的 解 ( 511 ) 对 :进行 微分 


DA EN nit 
dt dt dt 


=-Avsinvt + Bvcosvt 


2 
eet) -Av Biase + By deosrt 
dt? dt dt 





=-Arcosvt - Brsinvt 


代入 微分 方程 ( 5.32 ) 


{-Arcosvt - Br?sinvt} + 2y{ -Avsinvt + Bucosut} + w*(Acosut + Bsinzt) 
= fcosut 

进行 整理 
{(w?- 7)A + 2B} cosit + {2wA + (w*- 7)B}sinvt = foosut 

为 满足 上 式 


(u-)A+2wB=f 
-2wA+ (u?-A)B=0 


必须 成 立 。 对 这 个 并 列 方程 进行 求解 ， 则 常数 4 和 有 为 
sg wy 
(w=) + (2)? (P-P + (2? 
根据 假设 的 解 ( 5.33 ) 得 到 特 解 
oe WE ee A pig, 
(WP + (2)? (=P) + Qn)? 
最 后 ， 微 分 方程 (5.32 ) 的 通 解 为 齐 次 方程 48- + 21 SE + te =O 的 通 解 (5.16) 
Z(t)=e (acos Qt + bsin Qt), Q=4 u- F 
和 非 齐 次 方程 ( 5.32 ) 的 特 解 (5.34 ) 之 和 ， 因 此 
a(t) =e" (acos Qt + bsin Qt) 
+ els feos ut + 2w fsinvt 
(P-P + (2w)? (- 7+ (2? 


x(t) 





(5.34) 





(5.35) 
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Q xraswer. 








O 引入 外 力 的 情况 下 ， 挫 在 弹簧 上 小 球 运动 的 解 


图 有 外 力 情况 下 的 解 


© 考虑 小 球 在 原点 静止 的 状态 下 ， 有 周期 性 的 外 力作 用 于 小 球 ， 且 外 力 的 角速度 与 
固有 角速度 相等 。 也 就 是 说 Y=w。 初 始 条 件 {=0 时 ，z(0)=z。, [dz/dt],u=0。 首 先 ， 
通 解 (5.35 ) 中 代入 7=w 


a(t) = e*(acos Qt + bsin A) +—/ sinut (5.36) 
2w 
应 用 初始 条 件 z(0)=z。， 则 
z(0)=erafacos(2.0)+bsin(2.0)}+ Fiwon sani (537) 


另外 ， 将 解 (5.36) 对 + 进行 徽 分 
了 


至 = -ye "(acos Qt + bsin Qt) + e(-aQsin Qt + bQcos Qt) + pat 
ad 


= e"{(~7a +bA)cos Mt ~ (yb + aN)sin A} 二 ee 








因此 
加 = €77{(~y0 + 6Q)eos(20) - (yb+a)sin(20)} + Loostvo) 
t= 
=-7a+00+L00 (5.38) 
27 
最 后 根据 式 (5.37) 和 (5.38) 得 到 
a=0 a-0 
f “ f 
70+ bQ+L00 “on 
ee tay 22 (539) 
因此 ， 把 式 (5.39) 代入 通 解 (5.36) 中 ， 则 解 为 
ef -了 tg 
a(t) e pait) pg 
为 了 使 其 变 得 简单 ， 使 Y<<w， 则 Q ~ w， 解 可 以 表示 为 


$ 


a(t) = ——(1-e)sinut 
2w 


© w 
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有 周期 性 外 力 情况 下 的 解 振幅 逐渐 变 大 。 
(强迫 报 动 ) 











像 这 样 的 振动 ， NN | 5 


AN, REMEE 
与 扳 幅 恒定 的 阐 谐 报 


MRF, RH ARR 
物体 的 固有 角速度 相等 的 情况 
T, 振幅 达到 最 大 。 


这 种 现象 叫做 共振 ， 
也 叫 共 鸣 。 
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另外 ， 


















































阻力 小 的 情况 下 ， 
REE XK. 


车 有 一 定 程度 的 阻力 ， 
相 比 之 下 振幅 会 较 快 达 
到 恒定 值 。 








像 这 样 拉 伸 的 话 ， 
HARREN, 




















实际 上 ， 也 有 吊桥 VO 
在 风力 作用 下 例 填 £| QR, TURTHR, È 


的 例子 。 筑 物 和 滑 期 机 的 设计 中 也 
i 要 考虑 避免 共振 现象 引起 
的 破坏 。 


由 于 风力 作用 形成 了 
RH, READERS 
桥 体 发 生 了 共振 。 
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MRR 


BARE LAW BE 












wry ter er 


现 家 说 明 
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Q! 






aa] 


那么 ， 到 此 为 止 我 们 通 
aia ft 
RR. \ 












spiana, Be, Qa fut 


oY: (at Ble” 
ei) oe n+ Bel, Pa Pat 














/ 
D 
t 






L 《 

Ose Macon ts al + dant, \ a 
2-0-7 RY val 
iB 

19! 

LAI 

RT 

re 


X 








a(t) = Asin ut + Booswt 








那 是 因为 对 大 地 的 
也 许 是 因为 数学 之 神 您 A PRET, 
举 的 例子 易 于 理解 吧 。 


ROM GT RES BER 
知识 有 多 么 重要 ! 
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要 想 掌握 模型 化 的 关键 


之 处 ， 你 还 应 该 试 着 解 
REA AMT 


到 此 为 止 ， 你 的 愿望 
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以 后 要 加 油 啊 ! 








































































































































































































































































































































































有 ? 个 


ae et Fe 
Room 你 们 怎么 还 没 走 啊 ? 起 如 如 村 开间 
4 Meg AC 
meg. 









































神仙 姐姐 ， 不 好 好 了 
解 一 下 大 家 的 生活 怎 
么 行 呢 ? 
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好 ! 





| 这 次 ， 我 可 以 实现 一 下 


女神 的 愿望 了 ! 





CA TIENE, 










会 受到 上 天 的 惩 h 
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Hi 


1 咖啡 的 冷却 
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出. 咖啡 的 冷却 | 


现在 ， 人 们 能 轻易 地 喝 到 各 种 咖啡 了 ， 对 于 喜欢 咖啡 的 人 来 说 这 真是 件 愉快 的 事 。 
要 说 省 事 儿 的 话 ， 还 是 要 说 日 本 ， 现 在 全 日 本 到 处 都 有 的 钠 装 咖啡 自动 售 货 机 了 。 

可 是 ， 装 人 经 中 的 热 咖啡 会 慢 慢 地 冷却 。 可 也 有 从 自动 售 货 机 里 出 来 热 得 不 能 碰 的 
咖啡 ， 直 接 喝 的 话 ， 会 被 烫伤 。 那 么 ， 什 么 时 候 才 适 于 饮用 呢 ? 

热 咖啡 的 冷却 也 能 用 微分 方程 建 模 。 如 果 咯 啡 的 温度 比 气 高 的 话 ， 咖 啡 的 温 
度 就 会 下 降 。 研 究 一 下 温度 随时 间 的 变化 ， 就 可 以 明白 其 与 周围 空气 的 温差 成 正比 !。 
若 设 时 间 为 :， 咖 啡 的 温度 为 7 ( 1) ,周围 空气 的 温度 为 T， 比 例 常数 为 k， 则 


到 =- 区 全 描述 咖啡 的 温度 随时 间 变化 的 微分 方程 








这 就 叫做 牛顿 冷却 定律 2。 
我 们 对 这 个 方程 进行 求解 首先， 确认 一 下 变量 。 





两 边 分 别 进行 积分 





lex a+] a 


式 子 的 左右 两 边 分 别 被 变量 7 和 ;的 积分 成 功 分 离 。 


1 如 果 有 温度 计 、 腕 表 、 画 图 用 纸 的 话 ， 就 能 立即 确认 。 
2 该 法 则 说 明 的 不 是 温度 随时 间 的 变化 ， 而 是 热量 随时 间 的 变化 ， 如 果 热 容 不 随 温度 发 生变 化 的 话 ， 也 表示 同样 的 结果 。 





| ar-a | dt 
T-Te 


左右 两 边 分 别 求 积 分 ， 则 


Jar- ilr-Tl+c 


T-T. 
af dt=-kt+C 
综合 积分 常数 ， 则 
In|T- T=-kt+C 


温度 关于 时 间 的 函数 7 (+) ， 则 变 成 
TD-T=ew'c ”人 一 微分 方程 的 解 
到 此 ， 微 分 方程 被 解 出 来 了 。 
其 次 是 确定 积分 常数 。 假 设 初始 条 件 为 ，r=0 时 咖啡 温度 为 TO)=Tu(>7Jj， 根 据 式 
(2.31) 得 出 
TO)-Ti=ee 
则 
T,-T,+e 
那么 ， 描 述 咖啡 温度 的 函数 为 
T)=(T,-Te*+T, 

















RS, 现实 世界 | 数学 世界 TOST 


这 里 要 说 明 咯 啡 如 何 随时 间 流 逝 而 冷却 的 


© 咖啡 冷却 现象 的 说 明 
3 这 里 指 的 都 是 热 咖啡 。 
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咖啡 冷却 解 曲线 为 核 训 变 解 曲线 的 上 移 图 +。 当 参数 变 大 的 时 候 ， 咖 啡 会 加 速 冷 
却 。 


t 
加 


将 核 误 变 的 解 曲线 向 上 平移 就 能 得 到 嘿 啡 冷却 的 解 曲线 
© 表示 咖啡 的 温度 7 随时 间 / 变 化 的 曲线 


也 就 是 说 ， 过 热 的 咖啡 从 自动 售 货 机 出 来 的 时 候 ， 仅 仅 用 很 短 的 时 间 测 定 温度 的 变 
化 的 参数 k 的 话 ， 就 可 以 预测 到 什么 时 候 可 以 喝 了 。 但 是 ， 预 测 之 后 ， 就 不 能 呼 呼 吹 咖 
啡 了 。 因 为 上 发 生变 化 的 话 ， 预 测 也 会 发 生 偏离 哦 。 


4 如 果 较 为 熟练 的 话 ， 仅 仅 看 到 微分 方程 就 能 在 心里 想象 出 其 曲线 的 样子 。 只 有 加 强 练习 才能 做 到 。 


|2. 火箭 的 飞行 | 

有 一 个 位 于 太平 洋 岸 的 小 镇 叫 大 树 叮 ， 在 北海 道 的 十 胜地 区 ， 它 以 打造 “宇宙 之 
城 ” 为 目标 ， 进 行 混合 燃料 火箭 5 的 发 射 实验 ， 大 部 分 火箭 都 是 通过 推进 剂 与 氧化 剂 进行 
化 学 反应 ， 把 燃烧 的 气体 合 弃 在 后 边 飞 行 的 。 由 于 火箭 不 能 舍弃 所 有 自身 重量 ， 它 的 速 
度 是 有 上 限 的 。 这 个 速度 的 上 限 叫做 到 达 速 度 。 那 么 ， 我 们 试 着 计算 一 下 火箭 的 到 达 速 
度 。 

设 时 间 为 ,火箭 的 速度 为 V (+1) ,火箭 的 质量 为 M 从 火箭 的 马达 排 气 的 速度 为 v 的 
话 ， 微 分 方程 为 


uaa 描述 火箭 速度 的 微分 方程 


这 是 一 个 可 分 离 变 量 微分 方程 。 对 该 方程 进行 求解 ， 设 舍弃 质量 之 前 火箭 的 质量 为 M.， 
舍弃 质量 之 后 的 火箭 质量 为 M， 则 


Vevla( 7) 一 微分 方程 的 解 


这 个 方程 就 叫做 齐 奥 尔 科 夫 斯 基 方程 ， 利 用 它 能 够 计算 火箭 的 到 达 速 度 。 比 如 说 ， 
使 用 排 气 速度 为 v=1500m/s 的 火箭 马达 ， 发 射 之 前 火箭 的 重量 为 M=10kg， 推 进 材料 燃烧 
后 火箭 的 重量 为 M=9kg 的 情况 下 ， 小 型 火箭 的 到 达 速度 为 158m/s。 


5 在 推进 剂 与 氧化 剂 中 使 用 固体 与 液体 的 火箭 。 
6 就 像 离子 推进 火箭 那样 ， 使 用 电能 的 火箭 也 在 慢 慢 地 得 到 实际 的 应 用 。 
7 火箭 就 像 运载 重 物 的 卡车 一 样 ， 因 此 火箭 本 身 与 重 物 的 质量 不 可 能 什么 都 不 留 下 。 
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[ 3. 心理 量 | 
我 们 身体 反应 中 的 有 些 东西 也 可 以 用 可 分 离 变量 微分 方程 来 说 胃 。 比 如 说 ， 闭 上 腿 
贱 ， 拿 过 来 一 个 108 的 物体 。 当 质量 变 成 两 售 20g 的 时 候 ， 感 觉 不 出 两 人 的 重重， 感觉 到 
两 信 重 量 的 时 候 ， 其 实 全 的 应 该 是 1008 左 右 重 量 的 东西 了 。 我 们 的 感觉 也 是 亲 条 对 数 本 
数 变 化 的 。 不 仅仅 是 对 重量 的 感觉， 对 声音 大 小 的 感 沉 等 也 是 如 此 。 这 就 叫做 费 荐 纳 定 
it. 
设 物理 性 强度 为/， 感 性 强度 为 E0，1 为 常数 的 话 ， 微 分 方程 可 以 写成 
AELE -一 措 述 感觉 性 强度 的 微分 方程 
ad I 
进行 求解 ， 则 
工 
B=Kin(7) 。 一 微分 方程 的 解 


0 


这 里 ， 我 们 设 与 能 够 感受 到 的 最 低 刺 激 相对 应 的 反应 为 0。 
EL)=0 
参数 4 和 /根据 人 的 不 同感 觉 而 异 。 


8 也 有 丝 感 的 人 和 迟钝 的 人 。 


[4 广告 的 效果 | 


首先 求 一 下 较 容易 求 得 的 解 ， 然 后 对 其 进行 补充 ， 使 其 满足 原来 想 要 求解 的 微分 方 
程 ， 这 种 解法 叫做 常数 变易 法 ， 它 不 仅 对 于 自然 现象 ， 对 于 社会 现象 也 是 很 有 效 的 。 

世上 充满 了 五 花 八 门 的 广告 。 为 什么 会 有 如 此 多 的 广告 呢 ? 这 是 因为 ， 不 宣传 的 
话 ， 销 售 额 会 下 降 "。 实 际 上 ， 如 果 完全 没有 促销 的 话 ， 销 售 状况 会 呈 指 数 函数 式 地 减 
少 "。 我 们 把 这 种 情况 用 微分 方程 进行 模型 化 。 

设 入 和 /为 常数 ， 根 据 调查 结果 ， 销 售 速度 ( 单位 时 间 内 出 售 的 数量 ) 8 随时 间 ; 的 
变化 能 用 下 面 的 指数 函数 进行 说 明 。 

Sse" (A1) 

当然 ， 在 只 有 一 个 人 感 兴趣 的 情况 下 ， 只 有 在 概率 上 能 知道 这 个 人 到 底 买 还 是 不 
买 。 那 么 ， 在 很 多 人 感 兴趣 的 情况 下 ， 不 做 广告 时 就 可 以 用 式 ( A.1 ) 所 描述 的 函数 来 
说 明 销售 速度 。 


t 
o 


多 不 做 广告 的 情况 下 ， 销 售 速度 随时 间 的 变化 


这 个 函数 可 以 作为 可 分 离 变 量 微分 方程 的 解 ， 早 已 为 我 们 所 熟知 。 参 考 可 分 离 变 量 
微分 方程 的 情况 ， 满 足 这 个 函数 的 微分 方程 可 以 表示 成 ! 


这 就 是 不 做 广告 的 情况 下 销售 速度 的 微分 方程 。 


9 或 者 因为 害怕 销售 阁下 滑 。 
10《 用 微分 方程 做 数学 模型 》(Burghes, David; Borrie, Morag 著 ， 日 本 评论 社 )P76。 数 据 有 些 陈旧 。 
11 对 真 假 有 所 怀疑 的 人 ， 请 把 函数 代入 进行 确认 。 
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那么 ， 做 广告 的 话 会 变 成 什么 样 呢 。 做 广告 的 话 ， 看 见 广告 的 人 当中 会 出 现 购买 
商品 的 人 =。 看 见 广告 的 人 当中 有 恒定 比例 的 人 购买 商品 的 话 ， 销 售 速度 应 该 与 看 见 广 
告 的 人 的 比例 4 (1) 呈正 比 地 增加 。 但 不 能 希望 做 了 广告 ， 一 定 就 能 卖 很 多 。 由 于 能 够 
购买 商品 的 人 数 有 限 ， 不 久 便 会 达到 项 点， 销售 速度 也 会 变 成 恒定 值 8?。 这 也 就 是 说 ， 
销售 速度 有 极限 M ( 销售 速度 的 最 大 值 ) , 当 实际 销售 速度 $ 接 近 销 售 速度 的 最 大 值 M 的 
时 候 ， 销 售 速度 $ 会 变 迟 缓 。 也 就 是 说 ， 销 售 速度 $ 与 销售 速度 的 最 大 值 M 和 销售 速度 $ 
之 差 (销售 速度 上 升 的 空间 ) 呈正 比 。 速 度 最 大 值 M 对 应 的 销售 速度 上 升 空间 M-5 之 比 
(MS) /M 表 示 的 是 销售 速度 增长 的 比例 。 综 上 所 述 ， 设 WH, 那么， 做 广告 之 后 
销售 速度 随时 间 的 变化 与 
M-S 
M 
成 正比 ， 因 此 微分 方程 可 以 表现 为 


aS Age ya MS 
dt M 


右边 稍微 进行 整理 ， 则 


m-s AS A 
-AS+ 9A Ao a AS e774 a OAST) S494 


yA: 


则 微分 方程 


ds 4 
I (asri) 9+y4 ORERE (A2) 


就 变 成 了 非 齐 次 方程。 立刻 用 常数 变易 法 进行 求解 
与 非 齐 次 方程 (A.2 ) 对 应 的 齐 次 方程 为 

ds A 

对 ~-(X+7 咎 )5 一齐 次 方程 


首先 用 常数 变易 法 对 其 进行 求解 。 变 量 分 离 ， 则 


(A3) 


12 不 能 断定 只 要 做 了 广告 就 会 出 现 购买 的 人 ， 这 句 话 的 意思 是 说 假设 这 个 广告 是 一 个 有 效 的 广告 。 
13 卖 贷 的 人 想 扩大 市 场 ， 但 是 不 论 多 想 扩大 销售 规模 ， 还 是 被 地 球 上 的 人 口 数量 所 限制 。 
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分 别 对 两 边 进行 积分 
blsl--(A+y 总 ji+c 
所 求 的 解 就 变 成 


seattle 


把 常数 综合 为 C-+e 则 
Saco TH! -FENER (a4) 


为 齐 次 方程 ( A.3 ) 的 通 解 。 


其 次 ， 对 齐 次 方程 的 通 解 ( A.4 ) 进行 补充 使 其 满足 非 齐 次 方程 ( A.2 ) 。 把 常数 
替换 为 时 间 : 的 函数 ， 则 


sonal) BUENE RR 


由 于 这 个 假设 的 解 满足 微分 方程 ( A.2 ) ， 代 入 式 中 进行 微分 
A 


-(rA)e (rere -(ayA)e 
w, ( ti) set (-(a+7Ae ( it) J-ren ( nie) +yA 


HAEN 

an __ , bor) 

dt yAe 

替换 成 a= A+ vA TR, 则 

而 = 证 Í P RE 

a 

所 求 的 通 解 为 

5 的 -74 之 +0) ew roent FANRAN (As) 

a a 


我 们 确定 一 些 具体 的 条 件 来 看 一 下 这 个 解 是 如 何 说 明 广告 效果 的 。 我 们 认为 ， 即 使 
不 做 广告 也 会 有 一 定 比例 的 人 购买 商品 。 因 此 开始 做 广告 的 时 间 设 为 :=0， 把 销售 速度 5 
的 初始 值 设 为 $ (0) =5,。 另 外 ，0<+<7 的 期 间 内 ， 做 一 定 的 广告 ， 看 到 广告 的 人 的 比 
例 是 恒定 的 。 
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假设 
A(W=A(0<1<D) 

把 销售 速度 的 初始 值 5 (0) =5, 代 入 解 (A.5 ) , 则 
S(0)=5,= r +cerea 


那么 , 0<1<7 期 间 内 的 销售 速度 为 
so- a(s- A) 一 0<t<7 期 间 内 的 销售 速度 (A6) 


广告 刚 开始 出 现 的 时 候 ， 虽 然 销售 速度 会 急剧 增加 ， 但 是 速度 的 增加 会 逐渐 地 变 迟 
缓 ， 然 后 接近 恒定 的 值 。 


S, 
t 
o T 


多 做 广告 的 情况 下 ， 销 售 速度 随时 间 的 变化 


如 果 停止 做 广告 ， 会 变 成 什么 样 呢 ? t=7T 时 停止 广告 的 话 ，t>T 时 由 于 4=0， 因 此 
微分 方程 (A.5 ) 为 


通 解 与 调查 结果 所 得 式 子 ( A.1 ) 相同 ， 则 
Szeto 人 熏 齐 次 方程 的 通 解 (A7) 
但 是 由 于 做 广告 之 后 销售 速度 持续 增长 ， 这 次 就 从 t=7 时 刻 的 销售 速度 8 (7) 开始 
减少 。 把 9<t<7 期 间 的 销售 速度 (A.6) PRAET, W 


S(T)= (3, H) or 


因此 ， an (A7) , 则 


变 成 


S(T)=c'e t= A, (s,- £) er 
a a 


C= (4+ (5,- TA) at) ear 
最 后 ， 解 为 
S(t)= (4 R (5,- A) pot) owen 僵 t>T 时 的 销售 速度 
如 果 停止 广告 ， 销 售 速度 会 逐渐 地 减 小 。 
s 


M 
s) 


t 
了 | 


© 停止 广告 之 后 销售 速度 随时 间 的 变化 


虽然 现实 的 经 济 中 没有 那么 纯粹 的 事情 8&， 但 是 涉及 到 很 多 人 的 行为 的 时 候 ， 个 性 
被 平均 化 ， 人 们 的 行动 也 能 用 微分 方程 来 说 明 ， 真 是 好 有 趣 呀 。 


14 实际 上 ， 正 在 制作 更 为 复杂 的 模型 。 
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5. 积分 因子 解法 


第 3 章 中 ,我 们 看 过 了 用 常数 变易 法 求解 非 齐 次 线性 微分 方程 的 方法 ， 但 还 有 一 些 
其 他 解法 。 粗 略 地 看 一 下 
非 齐 次 方程 





到 plz)y= qlz) 。 全 非 齐 次 方程 
的 两 边 除 以 菜 个 函数 (a), M 


toy + pu(z)p(z)y = w(z)q(2) (A8) 
假设 ， 这 个 函数 & (zj 碰巧 有 


Sete) i 


的 性 质 ， 那 么 ， 这 个 微分 方程 的 右边 就 是 只 合 z 的 函数 ， 进 行 积分 ， 则 
way | wzjqzjdz (A9) 


这 样 ， 就 能 够 求解 了 。 
我 们 可 以 试 着 调查 一 下 ， 是 否 存在 那样 一 个 碰巧 的 函数 w(z)。 由 于 乘积 的 积分 为 


dc) _ dule) dy 
a a R 


与 微分 方程 (A.8 ) 相 比 


SHE) papa) 一 可 分 高 变量 微 分 方程 


如 果 上 述 关系 成 立 ， 那 就 变 得 非常 方便 了 。 这 是 可 分 离 变量 微分 方程 


求解 


In|u(@) f pz)dz 


“J(Z)=+er 二 积分 因子 
这 样 ， 所 得 到 的 这 个 很 方便 的 函数 叫做 积分 因子 。 了 解 了 积分 因子 ， 那 么 根据 式 
(A9), 则 


Joras 
H(z) 
也 就 是 说 ， 能 够 求 得 非 齐 次 方程 的 通 解 。 
处 理 广告 效果 的 非 齐 次 方程 (A2) 、 用 a= + Y -分 进 行 替换 
Sassa RIO 
能 够 试 着 利用 积分 因子 对 上 述 方程 进行 求解 。 积 分 因子 vb) 则 
H(t)=+e!' be C=+eCert=cent 
看 到 通 解 (A.10 ) 就 能 明白 ， 积 分 因子 也 有 常数 倍 的 任意 性 ， 包 含 积分 常数 的 
C= +e5 为 ] 也 没有 关系 。 也 就 是 说 
ute BOAT 
就 是 积分 因子 。 最 后 


Juterraae nau yA 
50= 一 一 一 = 一 一 一 -一 +cet +- 非 齐 次 方程 的 通 解 
u(t) et a 


就 变 成 通 解 。 正 好 与 通 解 (A.5 ) 一 致 。 若 发 现 用 积分 因子 求解 方便 的 函数 *， 那 就 能 够 
简单 地 求 出 非 齐 次 方程 的 解 。 


15 若 不 进行 积分 就 不 能 发 现 积分 因子 ， 因 此 不 是 什么 时 候 都 能 发 现 积分 因子 。 
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146 页 中 用 常数 变易 法 处 理 的 非 齐 次 方程 
Levande) FUNE AARE) 

只 要 稍微 进行 延展 的 话 ， 就 可 以 写成 这 种 形式 
yd 二 微分 方程 非 线性 微分 方程 ) (All) 


n=0 的 情况 为 非 齐 次 方程 、n=1 的 情况 为 齐 次 线性 微分 方程 。 除 此 之 外 的 情况 ， 就 
变 成 伯 努 利 微分 方程 的 非 线性 微分 方程 。 伯 努 利 微 分 方程 为 通过 实施 变量 变换 ， 得 到 线 
性 微分 方程 ， 能 够 用 常数 变易 法 进行 求解 。 





© 齐 次 方程 、 非 齐 次 方程 和 伯 努 利 方程 


伯 努 利 方程 通过 置换 因 变 量 y 
ze 变量 变换 (A12) 
六 
能 够 将 其 变 成 线性 微分 方程 。 我 们 试 一 试 。 伯 努 利 方程 ( A.11 ) 的 两 边 除 以 z", 则 
1 d 1 
六 (A.13) 


MA, 式 (A.12 ) 的 两 边 对 z 进 行 微分 ， 则 
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将 上 式 子 与 式 (A.12 ) RARR (A.13 ) 中 ， 则 可 以 写成 
-人 +o- 
为 了 更 容易 理解 ， 两 边 同 时 乘 以 1-n， 则 


= +(n-I)p(a)z=(n-1)g(z) 全 非 齐 次 微分 方程 (线性 微分 方程 ) 


就 变 成 线性 微分 方程 了 。 这 样 一 来 ， 就 方便 求解 了 。 
车 是非 齐 次 微分 方程 ， 用 常数 变易 法 和 积分 因子 法 都 能 够 求解 。 也 就 是 说 ， 伯 努 利 
微分 方程 ， 是 通过 变量 替换 、 进 行 线性 化 的 方法 来 求解 的 。 
第 107 页 中 用 分 解 部 分 分 数 这 种 分 离 变量 的 方法 所 解 的 物流 模型 的 微分 方程 ， 实 际 
上 就 是 伯 努 利 微分 方程 。 我 们 来 想 一 下 ， 人 口 用 P 表 示 ， 人 口 的 饱和 量 用 表示 ， 马 尔 
萨 斯 半径 用 /表示 ， 那 么 就 存在 这 个 微分 方程 
a l-g)” 


把 它 的 右边 展开 、 变 形 ， 得 到 


dP Pa p? 
Set. 僵 伯 努 利 微分 方程 ( 非 线性 微分 方程 ) (A.14) 


这 个 微分 方程 立刻 变 成 了 n=2 的 伯 努 利 微分 方程 。 
非 线性 微分 方程 ( A.14 ) 的 两 边 同时 除 以 PP， 得 到 


eA TE 

Pa "POE (A.15) 
这 里 ， 把 因 变量 P 做 变量 替换 

z4 < 变量 替换 (A.16) 
那么 ， 对 t 求 微分 就 得 到 

dz LP 

w Fra 
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利用 上 式 ， 非 线性 微分 方程 ( A.15 ) 


Beye ERON EMARE) (a17) 


就 变 成 线性 微分 方程 了 ， 这 就 是 非 齐 次 微分 方程 。 之 后 ， 用 常数 变易 法 进行 求解 。 首 
Fe, 求解 齐 次 微分 方程 


dz 
> 0 僵 齐 次 微分 方程 


变量 分 离 





僵 齐 次 方程 的 解 (A.18) 

a(t) = te= ce 

假设 C= + ec。 设 系数 c 为 时 间 ! 的 函数 

C=C() 
把 微分 方程 的 解 (A.18 ) 写成 

z(t)=c(t)e” 
将 其 代入 所 有 满足 的 非 齐 次 微分 方程 

d(c(t)e™*) 


tje” = 
dt +uc(t)e 


x 


dod) ey de oy was 
dt eet a c(t) +uc(t)e 


Seo. ev — eMe(t) tuc(the = 


然后 进行 积分 
得 到 
e(=4 j etdtec= +e 
K K 
因此 ， 解 为 
a=(S+e}em 一 非 齐 次 方程 的 角 
利用 式 (A.16 ) ,把 z 恢 复 为 P， 则 解 为 


M2 TE A ee 僵 就 是 微分 方程 的 解 (A19) 


与 前 一 章 中 一 样 ，t=0 时 的 人 口 为 P。 
Ke? K 


o PK 1+Ke’ 





把 上 式 代入 解 (A.19 ) 中 
Ke _ KP, 
K-P —Py)et+ P. 
e+K (a) (K-Pye o 
与 108 页 的 (3.7) 一 致 %。 
条 条 大 道 通 罗马 ， 求 解 微分 方程 也 可 以 用 各 种 各 样 的 方法 。 首 先 ， 只 有 掌握 了 成 熟 
的 技术 ， 才 能 进一步 发 展 。 伯 努 利 微分 方程 经 过 线性 化 之 后 ， 就 能 够 通过 常数 变易 法 进 
行 求解 。 首 先 ， 求解 一 个 较 容易 的 齐 次 微分 方程 ， 然 后 对 这 个 解 进行 补充 使 其 满足 本 来 
要 求解 的 非 齐 次 微分 方程 。 


P(t)= 





为 了 闪 到 目的 





16 从 P 项 的 齐 次 方程 的 解 中 能 够 通过 常数 变易 法 求解 ， 也 能 通过 积分 因数 进行 求解 。 
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